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Abstract 

New results for the complete one-loop contributions to the masses and mixing effects 
in the Higgs sector are obtained for the MSSM with complex parameters using the 
Feynman-diagrammatic approach. The full dependence on all relevant complex phases 
is taken into account, and all the imaginary parts appearing in the calculation are 
treated in a consistent way. The renormalization is discussed in detail, and a hybrid 
on-shell/DR scheme is adopted. We also derive the wave function normalization factors 
needed in processes with external Higgs bosons and discuss effective couplings incorpo- 
rating leading higher-order effects. The complete one-loop corrections, supplemented 
by the available two-loop corrections in the Feynman-diagrammatic approach for the 
MSSM with real parameters and a resummation of the leading (s)bottom corrections 
for complex parameters, are implemented into the public Fortran code FeynHiggs 2.5. 
In our numerical analysis the full results for the Higgs-boson masses and couplings are 
compared with various approximations, and CT^-violating effects in the mixing of the 
heavy Higgs bosons are analyzed in detail. We find sizable deviations in comparison 
with the approximations often made in the literature. 
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1 Introduction 



A striking prediction of models of supersymmetry (SUSY) [1] is a Higgs sector with at least 
one relatively light Higgs boson. In the Minimal Supersymmetric extension of the Standard 
Model (MSSM) two Higgs doublets are required, resulting in five physical Higgs bosons: the 
light and heavy CP-even h and H, the CP-odd A, and the charged Higgs bosons H^. The 
Higgs sector of the MSSM can be expressed at lowest order in terms of Mz (or Mw), Ma 
(or Mh±) and tan (3 = ^2/^1, the ratio of the two vacuum expectation values. All other 
masses and mixing angles can therefore be predicted. Higher-order contributions give large 
corrections to the tree-level relations. 

The limits obtained from the Higgs search at LEP (the final LEP results can be found in 
Refs. [2,3]), place important restrictions on the parameter space of the MSSM. The results 
obtained so far at Run II of the Tevatron [4-6] yield interesting constraints in particular in 
the region of small and large tan (the dependence on the other MSSM parameters has 
recently been analyzed in Ref. [7]). The Large Hadron Collider (LHC) has good prospects 
for the discovery of at least one Higgs boson over all the MSSM parameter space [8-10] (see 
e.g. Refs. [11,12] for recent reviews). At the International Linear Collider (ILC) eventually 
high-precision physics in the Higgs sector may become possible [13-15]. The interplay of the 
LHC and the ILC in the MSSM Higgs sector is discussed in Refs. [16, 17]. 

For the MSSM with real parameters (rMSSM) the status of higher-order corrections to 
the masses and mixing angles in the Higgs sector is quite advanced. The complete one-loop 
result within the rMSSM is known [18-21]. The by far dominant one-loop contribution is the 
0{at) term due to top and stop loops {at = /i^/(47r), ht being the top-quark Yukawa cou- 
pling). The computation of the two- loop corrections has meanwhile reached a stage where 
all the presumably dominant contributions are available [22-36]. In particular, the 0{atas), 
C(q!j), 0{ai)as), 0{atah) and 0{al) contributions to the self-energies are known for vanish- 
ing external momenta. For the (s)bottom corrections, which are mainly relevant for large 
values of tan/3, an all-order resummation of the tan /3-enhanced term of (9(a;,(a<j tan/3)"') 
is performed [37-40]. The remaining theoretical uncertainty on the lightest CP-even Higgs 
boson mass has been estimated to be below ~ 3 GeV [41-43]. The above calculations have 
been implemented into public codes. The program FeynHiggs [23, 44-46] is based on the 
results obtained in the Feynman-diagrammatic (FD) approach [22,23,34,41]. It includes all 
the above corrections. The code CPsuperH [47] is based on the renormalization group (RG) 
improved effective potential approach [26,35,36]. Most recently a full two-loop effective 
potential calculation^ (including even the momentum dependence for the leading pieces) has 
been published [49]. However, no computer code is publicly available. Besides the masses 
in the Higgs sector, also for the couplings of the rMSSM Higgs bosons to SM bosons and 
fermions detailed higher-order corrections are available [37-39,50,51]. 

In the case of the MSSM with complex parameters (cMSSM) the higher order corrections 
have yet been restricted, after the first more general investigations [52], to evaluations in the 

^In Ref. [48] the symmetry relations affecting higher-order corrections in the MSSM Higgs sector have been 
analyzed in detail. It has been shown for those two-loop corrections that are implemented in FeynHiggs 2 . 5 
that the counterterms arising from multiplicative renormalization preserve SUSY, so that the existing result 
is valid without the introduction of additional symmetry-restoring counterterms. It is not yet clear whether 
the same is true also for the subleading two- loop corrections obtained in Ref. [49]. 
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effective potential (EP) approach [53,54] (at one-loop, neglecting the momentum-dependent 
effects) and to the RG improved one-loop EP method [55,56]. The latter ones have been 
restricted to the corrections arising from the (s)fermion sector and some leading logarith- 
mic corrections from the gaugino sector^. Within the FD approach the one- loop leading 
corrections have been evaluated in Ref. [57] . Effects of imaginary parts of the one-loop contri- 
butions to Higgs boson masses and couplings have been considered in Refs. [58-60]. Further 
discussions on the effect of complex phases on Higgs boson masses and decays can be found 
in Refs. [61-64]. A detailed comparison between the two available computer codes for the 
cMSSM Higgs-boson sector, FeynHiggs and CPsuperH, will be performed in a forthcoming 
publication. 

In the present paper we present the complete one-loop evaluation of the Higgs-boson 
masses and mixings in the cMSSM (see Ref. [65] for preliminary results). The full phase 
dependence, the full momentum dependence and the effects of imaginary parts of the Higgs- 
boson self-energies are taken consistently into account. Our results are based on the FD 
approach using a hybrid renormalization scheme where the masses are renormalized on-shell, 
while the DR scheme is applied for tan/3 and the field renormalizations. The higher-order 
self-energy corrections are utilized to obtain wave function normalization factors for external 
Higgs bosons and to discuss effective couplings incorporating leading higher-order effects. We 
provide numerical examples for the lightest cMSSM Higgs boson, the mass difference of the 
heavier neutral Higgses and for the mixing between the three neutral Higgs bosons. We 
compare our results with various approximations often made in the literature. All results 
are incorporated into the public Fortran code FeynHiggs 2.5 [23,44-46]. 

The rest of the paper is organized as follows. In Sect. 2 we review all relevant sectors 
of the cMSSM. Besides the tree-level structure of the Higgs sector, the renormalization 
necessary for the one-loop calculations is explained in detail. In Sect. 3 the evaluation 
of the one-loop self-energies is presented. The determination of the Higgs-boson masses 
from the propagators and of wave function normalization factors and effective couplings is 
described. Our numerical analysis is given in Sect. 4. Information about the Fortran code 
FeynHiggs 2.5 is provided in Sect. 5, more details about installation and use are given in 
the Appendix. We conclude with Sect. 6. 

2 Calculational basis 

2.1 The scalar quark sector in the cMSSM 

The mass matrix of two squarks of the same flavor, and qj^, is given by 

fMl + ml + Ml cos 2/3( J| - Q,sl) m, X* 

'V ^<?^5 M|^ + m2 + M| cos 2/3Q,4 

with 

Xq = Ag- fi* {cot f3, tan (3}, (2) 

^The two-loop results of Ref. [49] can in principle also be taken over to the cMSSM. However, no explicit 
evaluation or computer code based on these results exists. 
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where {cot/3, tan/5} applies for up- and down-type squarks, respectively, the star denotes a 
complex conjugation, and tan/3 = v-z/vi. In Eq. (2) M|, are real soft SUSY-breaking 
parameters, while the soft SUSY-breaking trilinear coupling Ag and the higgsino mass pa- 
rameter can be complex. As a consequence, in the scalar quark sector of the cMSSM 
Nq + 1 phases are present, one for each Aq and one for /i, i.e. N^ + l new parameters appear. 
As an abbreviation we will use 

V'x, = arg [Xq) , i^A, = arg [Aq) . (3) 

One can trade ifAq for ipxq as independent parameter. 

The squark mass eigenstates are obtained by the unitary transformation 

1)^-^(1) w 

with 

The elements of the mixing matrix U can be calculated as 



I Ml + ml + Ml cos 2/3(J| - Qqsl) - m| 
= ^ p— (6) 



mqX* 



Ml + M| cos 2/3(/| - Qqsl) + - ml , I ml - ml 



<?2 V 91 92 



Here Cq = cos 6q is real, whereas Sq = e ^'^^i sin 6jj can be complex with the phase 

Vsi^ = - fx, = arg {X*q) . 
The mass eigenvalues are given by 
1 



(7) 



M| + M|^ + /|M| cos 2/3 (9) 



2 2 I 

m^, ^ = rn,, H — 

<?1,2 <? 2 

T y/[M| - + M| cos 2/3(J| - 2g,4)]2 + 4ml\Xq\^ ] , (10) 
and are independent of the phase of Xq. 

2.2 The chargino / neutralino sector of the cMSSM 

The mass eigenstates of the charginos can be determined from the matrix 

y-^f M2 V2smp Mw\ .... 

\V2cos(3Mw 1^ J' 

In addition to the higgsino mass parameter fi it contains the soft breaking term M2, which 
can also be complex in the cMSSM. The rotation to the chargino mass eigenstates is done by 
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transforming the original wino and higgsino fields with the help of two unitary 2x2 matrices 
U and V, 



xt 



= V 




u 



'12) 



These rotations lead to the diagonal mass matrix 



m-± 
m-± 

A2 



U*XV"f. 



(13) 



From this relation, it becomes clear that the chargino masses m -± and m -± can be determined 
as the (real and positive) singular values of X. The singular value decomposition of X also 
yields results for U and V. 

A similar procedure is used for the determination of the neutralino masses and mixing 
matrix, which can both be calculated from the mass matrix 



Y 



/ Ml 


-Mz cos/5 
V Mz sin (3 





M2 
Mz c„ cos /3 
Mz Cw sin (3 



-Mz cos (3 
Mz Cw cos /3 


-/i 



Mz sin 13 \ 
Mz sin P 

/ 



(14) 



This symmetric matrix contains the additional complex soft-breaking parameter Mi. The 
diagonalization of the matrix is achieved by a transformation starting from the original 
bino/wino/higgsino basis. 



















\ 


x'2 

xl 


= N 


m 


















\xV 








V 











N*YN 



(15) 



The unitary 4x4 matrix N and the physical neutralino masses again result from a numerical 
singular value decomposition of Y. The symmetry of Y permits the non-trivial condition 
of using only one matrix N for its diagonalization, in contrast to the chargino case shown 
above. 



2.3 The cMSSM Higgs potential 

The Higgs potential Vh contains the real soft breaking terms ml and (with mf = mi + |/ip, 
= fh^ + l/ip), the potentially complex soft breaking parameter m^g, and the U(l) and 
SU(2) coupling constants gi and (72: 

+ \{gl + gl){Hl^Hu - H;^H,,f + \gl\H*M\ (16) 
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The indices = {1, 2} refer to the respective Higgs doublet component (summation over 
i and j is understood), and e^^ = 1. The Higgs doublets are decomposed in the following 
way, 



Hi 

n2 



H22 



Vl + 



71^ 



^2 + 7l( 



(17) 



Besides the vacuum expectation values Vi and V2, eq. (17) introduces a possible new phase ^ 
between the two Higgs doublets. Using this decomposition, Vh can be rearranged in powers 
of the fields, 



V, 



H 



+ 1 



X1.X2) M, 



<t><f>xx 



(P2 

Xi 

\X2j 



+ 



'1 5 S^2 



+ 



where the coefficients of the linear terms are called tadpoles and those of the bilinear terms 



are the mass matrices M 



and Mrf 



The tadpole coefficients read 



= -\/2{m\vi - cos^VLl^z + \{9l + 9l){vl - vl)vi), 
T<p2 = -V2{mlv2 - cos^VLl^i - 1(^1 + 9l){vl - ^^2)^2), 



T^^ = \/2sin^'|m^2k2 = —T. 



X2 



V2 



(19a) 
(19b) 
(19c) 



12) 



with ^' = ^ + arg(m 

With the help of a Peccei-Quinn transformation [66] /i and can be redefined [67] 
such that the complex phase of vanishes. In the following we will therefore treat 
a real parameter, which yields 



12 1 



12? 



e' = e 



The real, symmetric 4x4- matrix M^^^^ 
the following elements. 



and the hermitian 2x2-matrix 



M, 



M 



(l>X 



<P<t>XX 



M 



t 

<t>x 



(20) 



contain 



(21a) 



M 



-cos^mj^ - U9l + 9i)viV2 





— sin ^ m 



sin ^ 



2 

12 



- cos ^mj^-l {gf + 9^)viV2'' 
ml + \i9! + 9l)i^vl-vl) 



(21b) 
(21c) 
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'mj + l{gf + g^){vj-v. 
— cos C, 



- cos ^ mf 2 
ml + Hgf + gi)ivi 



12 



m 



12 



|5'i^i^2 



(21d) 



(21e) 



The non- vanishing elements of M^^ lead to CP- violating mixing terms in the Higgs potential 
between the CP-even fields (pi and 02 and the CP-odd fields xi ^'^'i X2 if 7^ 0- The mass 
eigenstates in lowest order follow from unitary transformations on the original fields, 



H 

A 

\GJ 



U 



n{0) ■ 



Xi 

\X2/ 



U 



c(0) ■ 



(22) 



The matrices Un(o) and Uc(o) transform the neutral and charged Higgs fields, respectively, 
such that the resulting mass matrices 



M 



diag 
hHAG 



and M' 



diag 
H±G± 



c(0)iVl<^±<^±U^(Q) 



(23) 



are diagonal in the basis of the transformed fields. The new fields correspond to the three 
neutral Higgs bosons h, H and A, the charged pair and the Goldstone bosons G and 

The lowest-order mixing matrices can be determined from the eigenvectors of M,^,^^^ and 
M0±0±, calculated under the additional condition that the tadpole coefficients (19) must 
vanish in order that vi and V2 are indeed stationary points of the Higgs potential. This 
automatically requires ^ = 0, which in turn leads to a vanishing matrix M,^^ and a real, 
symmetric matrix M^i^i. Therefore, no CP- violation occurs in the Higgs potential at the 
lowest order, and the corresponding mixing matrices can be parametrized by real mixing 
angles as 



Un(0) 



/ — sin a cos a \ 

cos a sin a 

— sin /3n cos /3n 

y cos/3n sin/Sj^J 



U, 



c{0) 



— sin Pc cos Pc 
cos (3c sin (3c 



(24) 



The mixing angles a, (3^ and (3c can be determined from the requirement that this transfor- 
mation results in diagonal mass matrices for the physical fields. It is necessary, however, to 
determine the elements of the mass matrices without inserting the explicit form of the mixing 
angles and keeping the dependence on the complex phase ^, since these expressions will be 
needed for the renormalization of the Higgs potential and the calculation of the tadpole and 
mass counterterms at one-loop order. 



2.4 Higgs mass terms and tadpoles 

In order to specify our notation and the conventions used in this paper we write out the 
Higgs mass terms and tadpole terms in detail. The terms in Vh-, expressed in the mass 
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eigenstate basis, which are hnear or quadratic in the fields are denoted as follows, 



Vr 



H 



const. — Th ■ h — Th ■ H — Ta ■ A — Tq ■ G 



i {h, H, A, G) 



{H~,G~)- 



m 



( ml 




TTT'hA 


T^hG\ 






miff 




^HA 


^HG 




H 




^HA 


m\ 


^AG 




A 


\Kg 




t^\g 


ml J 








2 

f^H-G+ 



















(25) 



Our notation for the Higgs masses in this paper is such that lowest-order mass parameters 
are written in lower case, e.g. m\, while loop-corrected masses are written in upper case, e.g. 
Ml 

For the gauge- fixing, affecting terms involving Goldstone fields in Eq. (25), we have cho- 
sen the 't Hooft-Feynman gauge. In the renormalization we follow the usual approach where 
the gauge-fixing term receives no net contribution from the renormalization transformations. 
Accordingly, the counterterms derived below arise only from the Higgs potential and the ki- 
netic terms of the Higgs fields but not from the gauge-fixing term. 

In order to perform the renormalization procedure in a transparent way, we express the 
parameters in Vh in terms of physical parameters. In total, Vh contains eight independent 
real parameters: Vi, V2, gf, g^i mt\, m^, ^^^d ^, which can be replaced by the parameters 
Mz, Mw, e, mH± (or m^), tan/3, Th, Th and T4. Thereby, the coupling constants gi and 
g2 are replaced by the electromagnetic coupling constant e and the weak mixing angle 9^ in 
terms of = cos6'w = M^/Mz, = ^yl - c^, 



e = gic^ = g2 Sw, 
while the Z boson mass Mz and tan/? substitute for vi and V2'- 

Ml = \{gl + gl){vl + vl), tan/3 = ^ 
The W boson mass is then given by 



M'w 



2/ 2 I 2\ 

2^2K +^'2) 



The tadpole coefficients in the mass-eigenstate basis follow from the original ones 
transforming the fields according to Eq. (22), 

Th = V2[—m\vi cos a —m\v2sma + cos ^ ("^i sin a +f2Cosa) 

- \{9l + gDi.'vl-vl) iy I cos a -t;2sina)), 
Th = v^l+mlvi sin a — m\v2 cos a + cos ^ (^1 cos a — V2 sin a ) 
+ + 9l){vl - vl){vi sin a + wgcosa )), 
= -a/2 sin ^ m^2 (""i cos /?n + f 2 sin /3n) , 



(26) 
(27) 

(28) 
(19) by 

(29a) 
(29b) 
(29c) 
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G 



tan(/3 - (3r,)TA. 



(29d) 



Using Eqs. (26) - (29) the original parameters can be expressed in terms of e, tan /5, Mz, 
Mw, Th, Th, Ta and either the mass of the neutral A boson, m^i, or the mass of the charged 
Higgs boson, mH± (it should be noted that Eqs. (29a)-(29d) yield only three independent 
relations because of the linear dependence of Tq on Ta)- The masses and mH± are 
related to the original parameters by 



m\ = ml sin^ /3n + ml cos^ /5n + sin 2/9n cos ^ 

- cos2f3n\{9f + 9l){vf - I'D, 

mlj± = ml sin^ Pc + f^l cos^ Pc + sin 2/5c cos ^ m^2 

- cos 2/34 (5-? + 5'2)(^i - ^2) + ifl'2K cos /9c + V2 sin /3c 

Choosing as the independent parameter yields the following relations, 
cos/? 



^^1 

^2 

91 

92 



v/2 sin/5 

e 

e/cw 

-^M| cos(2/?) + sin^ p/ {cos\p - /?n)) 



+ 



eTh cos Pn 
2c^s^Mz 

cTh cos /3n 



-M| cos(2/3) + m\ cos^ /?/ (cos2(/? - /?n)) 



2cwSwMz 
cos(2/3) - 

eTff sin P^ 

. 2c^s^Mz 
eTh sin /5n 



cos /5 cos Pn sin a + sin /5(cos a cos /?„ + 2 sin a sin /?„)) 
cos(a + P) cos /3n + 2 cos a sin P sin /5n) / (cos^(/? — /?n)) 



+ 



"^?2 

sin^ 
cos^ 

where 



fn 



2cwSwMz 

fm/ \/ fm ^ fs: 



-m^ sin 2/3 



sin a sin P sin /3n + cos/3(2 cos /3n sin a — cos a sin /?n)) 
(2cosQ;cos/?cos/3n + sin(a + /3) sin/3n) / (cos^(/3 — /3n)) 



[cos(/3 + a) + cos(/? — a) cos(2/3n)) 



+ 



eT, 



4c„SwMz 
cTa 



Ac^s^Mz 

sin(/3 + a) - sin(/3 - a) cos(2/3n)) / (cos^(/3 - /3n)) 



(30a) 
(30b) 

(31) 

(32) 

(33) 
(34) 



(35) 



(36) 

(37) 
(38) 
(39) 



2swCwMzcos(/3 - Pn) 



(40) 
(41) 
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We now give the bilinear terms of Eq. (25) in this basis, expressed in terms of m^i or 
mH±, depending on which parameter leads to more compact expressions. For the charged 
Higgs sector this yields, apart from mfj± itself, 

m\^Q+ = —m\± tan(/3 — (3^) (42a) 
g 

~ ^7i7 sin(a - /3c)/ cos(/3 - /3c) 

g 

~ 7777 cos(a - /3c)/ cos(/3 - /3c) 

g 

2MzSwCw 

^G-H+ = {'^l-G+y^ (42b) 
mg± = m^± tan^(/3 — /3c) 

-THCOs(a + /3-2/3c)/cos^(/3-/3c) (42c) 



2Mzs^c^ 
e 



+ TTTT^ sin(a + /3 - 2/3c)/ cos2(/3 - /3c). (42d) 



The neutral mass matrix is more easily parametrized by m^, as can be seen from the 
2x2 sub-matrix of the A and G bosons: 

= —m\ tan(/3 — /3n) (43a) 
g 

~ 77i7 sin(a - /3n)/ cos(/3 - (3^) 

g 

~ ^ ~ /^n)/ COs(/3 - /3n), 

m% = m\ tan^ (/3 - /3n) (43b) 
Tj^ cos(a + /3 - 2/3„)/ cos2(/3 - /3n) 



2MzSwC 



+ 7717^ sin(a + /3 - 2/3„)/ cos2(/3 - /3n). (43c) 

The CP- violating mixing terms connecting the h-/H- and the yl-/G-sector are 

^Ia = 7777 ^^^(" ~ ^'i)/ ~ ^°)' (44a) 

^hG = 1777 cos(a - /3n)/ cos(/3 - /3n), (44b) 

= -"^L. (44c) 

= 1777^ Ta sm{a - /3n)/ cos(/3 - /3n). (44d) 

Finally, the terms involving the CP-even h and H bosons read: 

ml = Ml sin^(a + (3) (45a) 
-|- m\ cos^(a — /3)/ cos^(/3 — /3n) 

+ 7777 '^os(" ~ sin^(" - /^n)/ cos^(/3 - /3n) 
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e 1 

-Th- sin(a — /3n)(cos(2a; — /9 — P^) + 3cos(/5 — Pn))/ cos^{P — /5n), 



m^j^ = -M| sin(a + P) cos(a + /5) (45b) 
+ sin(a — /5) cos(a — P)/ cos^{P — Pn) 

-TH-sin(a — /3) sin^(a — P^j cos^(/9 — /3n) 



-T/iCos(a — /5) cos^(a — P^j cos^(/5 — /3n), 



W 

= M| cos^(a + /3) (45c) 
+ sin^(a — P)/ cos^(/9 — /?„) 
e 1 

TH'-cos(a — /3n)(cos(2tt — — /5n) — 3cos(/? — /5n))/ cos^(/5 — /?n) 



2MzSwCw 2 
~ 7777 sin(a - p) cos^(a - /3n)/ cos^(/3 - /9n)- 

2.5 Masses and mixing angles in lowest order 

The masses and mixing angles in lowest order follow from the minimization of the Higgs 
potential. As mentioned above, this leads to the requirement that the tadpole coefficients 
T{h,H,A} and all non-diagonal entries of the mass matrices in Eqs. (42)-(45) must vanish (the 
tadpole coefficient Tq vanishes automatically if = holds). In particular, the condition 

= implies that the complex phase ^ has to vanish, see Eqs. (38)-(41), so that the Higgs 
sector in lowest order is CP-conserving. As a consequence, the well-known lowest-order 
results of the rMSSM are recovered from Eqs. (42)-(45). 

It follows from Eqs. (42a) and (43a) that the mixing angles have to obey 

Pc = Pn = p. (46) 

The lowest-order results for the Higgs masses can in principle be obtained from Eqs. (45a) 
and (45c) after the mixing angle a has been determined from Eq. (45b) by requiring that 
the right-hand side of the equation vanishes. More conveniently the Higgs masses can be 
determined by diagonalizing the 2x2 matrix in the 0i-02 basis, which corresponds to the 
entries (45) of the matrix of the neutral Higgs bosons in Eq. (25) with a set to zero. The 
lowest-order masses read 



[mlml] = ^(^ml + MlT ^J {m\ + M|)2 - 4m^M| cos^ 2/?j . 



(47) 



For the mixing angle a one obtains 

-{m\ + M|)sin/3cos/3 



a = arctan 



_M| cos^ p + m\ sin^ P - m\ 



^ < a < . (4J 



Finally, combining eqs. (30) and (46) relates the remaining masses tua and m^^i with each 
other, 

=m\ + clMl = m\ + M^. (49) 
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Depending on which of the masses and rriA is chosen as independent input parameter, 
the other mass can be determined from Eq. (49). Since the CP- violating mixing in the 
neutral Higgs sector implies that the CP-odd A boson is no longer a mass eigenstate in 
higher orders, we use the charged Higgs mass as input parameter for our analysis of 
the cMSSM. 



2.6 Renormalization of the Higgs potential 

We focus here on the renormalization needed for evaluating the complete one-loop correc- 
tions to the Higgs-boson masses and effective couplings (the latter corresponding to effective 
mixing angles) in the cMSSM. In our numerical analysis below we will also include two- 
loop corrections obtained within the FD approach, which up to now are only known for the 
rMSSM [22,23,41] (the renormalization of the relevant one-loop contributions is described 
in Refs. [22,23,41,42]). Also included will be the leading resummed corrections from the 
(s)bottom sector [37-39], which have been obtained in the cMSSM. 

In order to derive the counterterms entering the one-loop corrections to the Higgs-boson 
masses and effective couplings we renormalize the parameters appearing in the linear and 
bilinear terms of the Higgs potential. 



ML + SM, 



Th 
Th 
Ta 
tan P 



Th + 6Th, 
Th + 6Th, 
Ta + 6Ta, 
tan /? (1 + 5tan/3 ] 



(50) 



We express the counterterms in the mass-eigenstate basis of the lowest-order Higgs fields. 
While the parameter /3 is renormalized, the mixing angles (3^ and (3c (and also a) need not be 
renormalized. In carrying out the renormalization transformations it is therefore necessary 
to distinguish (3 from /5n and /3c (as we have done in (42)-(45)), i.e. Eq. (46) should only be 
applied after the renormalization transformations. 

For the counterterms arising from the mass matrices we use the definitions 



SMhHAG = Un(0) SM^^^^IJ 



'-<I>'PXX^ n{0) 



/ 6ml 




5^hA 


5Kg\ 






6mjj^ 


6m]jc 




S^HA 


6m\ 


6m\^ 






6m\(. 


6ml ) 






. 






^ 6m 







(51) 



(52) 



It should be noted that we need only seven independent counterterms in the Higgs sector, 
5m^±, 5M|, (5M^, 6Th, 6Th, 6Ta and 6 tan (3. This is due to the fact that in the expressions 
for the mass counterterms the renormalization of the electric charge, 6Ze, drops out at the 
one-loop level. Inserting the counterterms introduced in Eq. (50) and applying the zeroth 
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order relations T{h,H,A} = and /5n = /5c = in the coefficients of the first-order expressions 
yields for the CP-even part of the Higgs sector 

6ml = cos^(a - /?) + 5M| sm'^{a + (3) (53a) 

H {6Th cos(a — P) sin^(a — [3) + 6Th sin(a — /5)(1 + cos^(a — P))) 

2Mzs„c^ 

+ 5tan/3 sin (3 cos (3 {rriy^ sin 2 (a — /3) + M|sin2(a + /5)), 
StjiIh = l {Sml sin 2 (a - /3) - 5M| sin 2{a + p)) (53b) 

+ ^717^ ^^^^ - - ^^'^ - /^)) 

2M2S„Cw 

— (5tan/3 sin (3 cos /? (m^ cos 2(a — /?) + M| cos 2(a + /3)), 

= bm\ sin2(a - /3) + 5M| cos2(a + /?) (53c) 

{bTu cos(a — /3)(1 -|- sin^(a — /5)) + bT^ sin(a — /5) cos^(a — /?)) 

2MzSwCw 

— 5tan/3 sin (3 cos /? (m^ sin 2(a — /3) + M| sin 2(a + /?)), 

which has the same form as for the rMSSM. 
For the CP-odd part we obtain 

bni^j^Q = — — {—5Th sm{a — (3) — 5Th cos(a — (3)) — Staxifi m\ sin (3 cos (3 , (53d) 

'^"^G = 7777 (-(^^-ff cos(a - I3) + STh sin(a - /?)), (53e) 

ZMzS^C^ 

which again recovers the result of the rMSSM. 

For the counterterms arising from the CP-violating mixing terms we obtain 

^"T'hA = +7777 sin(a - p), (53f) 

^^Ig = +7777 cos(a - /?), (53g) 

^^HA = -^'^Ig^ (53h) 
5m]jG = 6mlA. (53i) 

Finally, the counterterms arising from the mass matrix of the charged Higgs bosons read 

6m]j-(.+ = — (STh sin(a - p) - 6Th cos(a ~ p) ~ i STa), (53j) 

zivizs^c^ 

— Stamp m'\j± sin P cos P , 

^^G~H+ = i^'mH-G+y^ (53k) 

5m^± = — — {—5Th cos(a — /5) + (5T/i sin(a — /?)). (531) 

2M^SwCw 

As mentioned above, we use mij± as independent input parameter. The counterterm 
SruA ill the formulas above is therefore a dependent quantity, which has to be expressed in 
terms of using 

5m\ = 5m]j± - 6M^, (54) 
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which follows from Eq. (49). 



For the field renormalization, which is necessary in order to obtain finite Higgs self- 
energies for arbitrary values of the external momentum, we choose to give each Higgs doublet 
one renormalization constant, 



In the mass eigenstate basis, the field renormalization matrices read 



(55) 



(h\ 

H 

A 
\GJ 



and 



/I + l6Zf,h \5ZuH 

\^ZhH 1 + \^Zhh 

\5ZhA \5Zha 

\ \5ZhG \5Zhg 



1 + \5Zh+h- 
\5Zg-h+ 

1 + \5Zh+h- 



bZ, 



hA 



1 + ^5Zaa 
\SZag 1 



\5ZhG 


\ 




\5Zhg 




H 


\5Zag 




A 


+ \5Zgg) 


\g) 



(56a) 















[g-I 






\5Zn-i 



Uz, 



G+G~ 







)■ 








)■ 


{G-j 



(56b) 
(56c) 



The renormalization according to Eq. (55) yields the following expressions for the field 
renormalization constants in Eq. (56): 



SZfj- 



H-G^ 



SZhh 


= sin^a + cos^a SZ-^^, 


(57a) 


5Zaa 


= sin^P 5Z-n-^ + cos^P SZ-^^, 


(57b) 


SZhH 


= sin a cos a {SZ-j-i^ — SZ-^J, 


(57c) 


SZag 


= sin (3 cos/5 {SZ-fi^ — SZy^J, 


(57d) 


SZhh 


= cos^a 5Z'Hi + sin^a SZ-j-i^, 


(57e) 


SZgg 


= cos^P 5Z-yi^ + sin^/5 5Z-yi^, 


(57f) 


SZh-h+ 


= sin^/5 6Zii^ + cos^/5 5Z-n^, 


(57g) 


SZg-h+ 


= sin P cos P {SZ-j-i^ — SZ-hJ, 


(57h) 


SZg-g+ 


= cos^P 5Z'n-^ + sin^P 5Zn^ . 


(57i) 



For the field renormalization constants of the CP-violating self-energies it follows, 

5ZhA = 5ZhG = 5Zha = 5Zhg = 0, 



(5^ 



which is related to the fact that the Higgs potential is CP-conserving in lowest order and 
Goldstone bosons decouple. 
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2.7 Renormalization conditions 



We determine the one-loop counterterms by requiring the following renormalization condi- 
tions. The SM gauge bosons and the charged Higgs boson are renormalized on-shell, 

ReSzz(M|) = 0, ReEw(M^) = 0, Re (M^±) = , (59) 

where the gauge-boson self-energies are to be understood as the transverse parts of the full 
self-energies. For the mass counterterms, Eq. (59) yields 

5M| = ReS^z(M|), ^M^, = ReT.ww{M^), 5m^± = ReSH+H-(M^±) . (60) 

It should be noted that Eqs. (59), (60) are strict one- loop conditions. Beyond one-loop order 
we define the mass of an unstable particle according to the real part of its complex pole, see 
the discussion in Sect. 3.4 below. 

The results for the self-energies can be decomposed as usual in terms of standard scalar 
one-loop integrals. Because of the appearance of complex phases, the coefficients of these loop 
integrals could in principle be complex. We have explicitly verified that this is not the case, 
i.e. the complex parameters appear in the results for the self-energies only in combinations 
where the imaginary parts cancel out. As a consequence, the only source for imaginary parts 
in the results for the self-energies are the loop integrals, as in the case of the rMSSM. 

As the tadpole coefficients are required to vanish, their counterterms follow from 

T{h,H,A}{l) + ST{h^H,A} = , (61) 

where T^h,H,A}{i) denote the one-loop contributions to the respective Higgs tadpole graphs: 

5Th = 6Th = -Th{i), 5Ta = -Ta(i)- (62) 

Concerning the field renormalization and the renormalization of tan/3, we adopt the 
DR scheme, 

6Zn, =6Z^ = - [Re |„=o] , (63a) 

6Zn, =6Z^ = - [Re S',, ^^^,] , (63b) 

(5tan/3 = ^ {SZn, - SZn,) = 6tanp^ (63c) 

i.e. the renormalization constants in Eqs. (63) contribute only via divergent parts. In 
Eqs. (63) the short-hand notation f'ijP') = d f {jP') / {d p^) has been used. As default value of 
the renormalization scale we have chosen in this paper yU,^^ = 

The DR renormalization of the parameter tan /5, which is manifestly process- independent, 
is convenient since there is no obvious relation of this parameter to a specific physical observ- 
able that would favor a particular on-shell definition. Furthermore, the DR renormalization 
of tan/3 has been shown to yield stable numerical results [19,45,68]. This scheme is also 
gauge-independent at the one-loop level within the class of gauges [68]. 

The field renormalization constants completely drop out in the determination of the 
Higgs-boson masses at one-loop order. They only enter via residual higher-order effects as 
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a consequence of the iterative numerical determination of the propagator poles described in 
Sect. 3.4 below. The DR scheme for the field renormalization constants is convenient in order 
to avoid the possible occurrence of unphysical threshold effects. Higgs bosons appearing as 
external particles in a physical process of course have to obey proper on-shell conditions. 
This issue will be discussed in Sect. 3.5. 

3 Higgs boson masses and mixings at higher orders 
3.1 Calculation of the renormalized self-energies 

At the one-loop level, the renormalized self-energies, S(p^), can now be expressed through 
the unrenormalized self-energies, S(p^), the field renormalization constants and the mass 
counterterms. As explained above, the counterterms arise from the Higgs potential and the 
kinetic terms, while the gauge-fixing term does not yield a counterterm contribution. The 
renormalization prescription of the gauge-fixing term induces counterterm contributions in 
the ghost sector, see e.g. Ref. [69] for further details. The counterterms from the ghost 
sector, however, contribute to the Higgs-boson self-energies only from two-loop order on. 



The renormalized self-energies read for the CP-even part, 

^hh{p^) = Sm(p^) + 5Zhh{p'^ - ml) - 5ml, (64a) 

thHip^) = S^h(p') + SZj^nip^ - \{ml + m^)) - 5m^^, (64b) 

^hh{p^) = T.hh{p^) + 5Zhh{p^ - mjj) - Smjj, (64c) 

and the CP-odd part, 

tAAip^) = ^AA^p") + ^Zaa^p" - m\) - 5m\, (64d) 

S^g(p') = S^g(p') + ^ZAcip" - Wa) - ^t^\g: (64e) 

Sgg(p') = Scg(p') + ^Zgcp" - 5ml. (64f) 

The CP- violating self-energies read (using Eq. (58)) 

tuAip") = ^hA{p') - Sml^, (64g) 

tf^cip') = S,g(/) - Smla, (64h) 

±HAip^) = Sha(/) - Sm^A, (64i) 

f:HG{p') = ^HG{p')-SmlG (64j) 

while for the self-energies in the charged sector one obtains 

t^H-H+{p^) = T.H-H+ip^) + 5Zh~h+{p^ - m]j±) - (5m^±, (64k) 

S//-G+(/) = S^-G+(p') + 5^H-G+(p' - >^±) - 5^l-G+^ (641) 

Sg-h+(/) = Sh-g+(p'), (64m) 

^G-G+{P^) = Sg-g+(p^) + SZg-g+p'^ - Sml±. (64n) 
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Figure 1: Generic Feynman diagrams for the /i, H , A, G self-energies (/ = {e, /i, r, s, 
6, M, c, t} ). Corresponding diagrams for the Z boson self-energy are obtained by replacing 
the external Higgs boson by a Z boson. 
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Figure 2: Generic Feynman diagrams for the iJ^, self-energies (/ = {e, fi, r}, d = {d, 
s, b}, u = {u, c, t} ). Corresponding diagrams for the W boson self-energy are obtained by 
replacing the external Higgs boson hj aW boson. 
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/ ^{e.M,r} /i,/2 H±,G^ xt,xt 




Figure 3: Generic Feynman diagrams for the h, H, A tadpoles (/ = {e, /i, r, d, s, b, u, c, 

0). 



The generic Feynman diagrams for the one-loop contribution to the Higgs and gauge- 
boson self-energies are shown in Figs. 1, 2. The one- loop tadpole diagrams entering via the 
renormalization are generically depicted in Fig. 3. As usual, all the internal particles in 
the one-loop diagrams are tree-level states. This implies in particular that diagrams with 
internal Higgs bosons do not involve CP-violating phases. The diagrams and corresponding 
amplitudes have been obtained with the program FeynArts [70] and further evaluated with 
FormCalc [71]. As regularization scheme we have used differential regularization [72], which 
has been shown to be equivalent to dimensional reduction [73] at the one- loop level [71]. 
Thus the employed regularization preserves SUSY [48,74]. 

In order to obtain accurate predictions for the Higgs-boson masses and mixings, in our 
numerical analysis below we will supplement the results for the one-loop Higgs self-energies 
in the cMSSM obtained in this paper with two-loop contributions where the dependence on 
the complex phases is partially neglected. The corresponding contributions will be described 
in Sect. 3.3. 

3.2 Special case: corrections to the charged Higgs-boson mass in 
the MSSM without CP-violation 

As a consequence of the mixing between the three neutral Higgs bosons in the presence of 
CP- violating phases in the Higgs sector, it is convenient to choose the mass of the charged 
Higgs boson, mH± , as the second free input parameter in the Higgs sector besides tan In 
the special case where the complex phases are zero (i.e. the rMSSM), on the other hand, 
one conventionally chooses the mass of the CP-odd Higgs boson, m^, as independent input 
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parameter instead of m/^±, so that the predictions for the neutral Higgs-boson masses do 
not involve charged Higgs-boson self-energies. 

In this case the mass of the charged Higgs-boson can be predicted in terms of the other 
parameters and receives a shift from the higher-order contributions. The results obtained 
in this paper can easily be applied to the special case of predicting the mass of the charged 
Higgs boson in the rMSSM, since the necessary ingredients are a subset of those entering 
the prediction for the neutral Higgs-boson masses in the cMSSM. 

The charged Higgs boson pole mass is obtained by solving the equation 

= 0, (65) 

where mf^± denotes the tree-level mass of the charged Higgs boson, Eq. (49), and T.n+H-ip'^) 
is defined in Eq. (64k). The mass counterterm, 5m?^±, is given in Eq. (54). In this approach, 
where is a free input parameter (m^ = Ma in our notation, since the tree-level mass mA 
does not receive higher-order corrections), the counterterm 6m\ can be fixed by the on-shell 
condition 

ReSAA(Ml) = 0, (66) 

leading to 

5m\ = ReT.AA{Ml) . (67) 

For earlier evaluations of the charged Higgs-boson mass, see Refs. [75,76]. A full one- loop 
calculation including a detailed numerical analysis can be found in Ref. [77]. 

3.3 Inclusion of higher-order corrections 

The numerical results for the Higgs-sector observables discussed below are based on the 
complete one-loop results obtained in this paper within the cMSSM, i.e. for arbitrary com- 
plex phases, supplemented by higher-order contributions. The renormalized self-energies are 
decomposed as 

S(p2) = S«(p') + S('n/) + -- - , (68) 

where S*^*^ denotes the contribution at the iih. order. 

In addition to the full one-loop contributions to S(p^), i.e. ^'•^-'(p^), in the cMSSM we in- 
corporate an all-order resummation of the tan /5-enhanced term of 0{ab{cxs tan (3)"^) including 
its phase dependence [37,38] . Since in the FD approach a result for the two- loop corrections in 
the t/t sector including the full phase dependence is not yet available (see however Ref. [61]) 
we take over the the leading two-loop QCD and electroweak Yukawa corrections obtained in 
the rMSSM [23,41], neglecting the explicit phase dependence at the two-loop level. All the 
contributions have been incorporated into the Fortran code FeynHiggs2.5 [23,44-46], see 
Sect. 5 below. 

3.4 Determination of the masses from the Higgs propagators 

In order to obtain the prediction for the Higgs masses beyond lowest order, the poles of the 
Higgs propagators have to be determined. Since the propagator poles are located in the 
complex plane, we define the physical mass of each particle according to the real part of the 
complex pole. 
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In determining the propagator poles one needs to take into account that the Higgs bosons 
mix among themselves, with the Goldstone bosons and with the gauge bosons. For the neu- 
tral Higgs bosons of the MSSM in the case with CV violation the Higgs propagators will in 
general receive contributions from the Higgs states /i, if, A, the Goldstone boson G, and the 
(longitudinal part of the) Z boson. The contributions of G and Z to the Higgs propagators 

appear from two-loop order on via terms of the form ^S^c^}^) j and [tl^zip^)^ , where 

(f) = h, H, A and S^^(p'^) = p^Tj^ziv^)- The contributions of G and Z are related to each 
other by the usual Slavnov-Taylor identities, ensuring a cancellation of the unphysical contri- 
butions. The mixing contributions with G and Z yield a sub-leading two-loop contribution 
(this contribution can be compensated at the propagator poles by a proper choice of the 
field renormalization constants, see e.g. Ref. [78]). As explained above, we supplement the 
one-loop Higgs-boson self-energies with the leading two-loop QCD and electroweak Yukawa 
corrections. Accordingly, the Higgs propagator terms induced by the mixing with G and 
Z are of the same order as terms that we neglect at the two-loop level. We will therefore 
neglect the effects induced by Higgs-boson mixing with G and Z in the determination of the 
Higgs-boson masses.^ Analogously, in the charged Higgs sector we neglect the mixing of 
with and W^. While the Higgs mixing with the Goldstone bosons and the gauge bosons 
yields subleading two-loop contributions to the Higgs-boson masses, it should be noted that 
mixing contributions of this kind can enter in Higgs decays or production processes already 
at the one- loop level (for more details, see Ref. [79]). 

According to the discussion above we can write the propagator matrix of the neutral 
Higgs bosons /i, if, A as a 3 x 3 matrix, ^hHAip^)- (The program FeynHiggs 2 . 5 allows to 
employ also the full 4x4 propagator matrix of all four scalar states H, A, G.) The 3x3 
propagator matrix is related to the 3x3 matrix of the irreducible vertex functions by 



^hHA{p') = - [ThHA{p')) , (69) 

where 

f/^//A(p')=^[p']l-M„(/)], (70) 

Mn(p^)= -thHip^) ml-±HH{p') -±ha{p') • (71) 

V -t,A{p') -^HAip') m\-tAA{p^)l 

Inversion of ThuAip^) yields for the diagonal Higgs propagators {i = h, H, A) 

An(p') = , (72) 

p'-m^ + ^fip') 

where Ahh^p"^), AnHip"^), Aaa{p^) are the (11), (22), (33) elements of the 3x3 matrix 
AhHAip"^), respectively. The structure of Eq. (72) is formally the same as for the case without 



■^We have explicitly verified that the numerical contributions of the mixing self-energies of the Higgs 
bosons with G and Z are indeed insignificant. 
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mixing, but the usual self-energy is replaced by the effective quantity S^f (p^) which contains 
mixing contributions of the three Higgs bosons. It reads (no summation over i,j, k) 

where the Tij^p"^) are the elements of the 3x3 matrix ThHAip"^) as specified in Eq. (70). 

For completeness, we also state the expression for the off-diagonal Higgs propagators. It 
reads {i 7^ j, no summation over i,j, k) 



A,,(p2) = ^^-^^ . . . (74) 

F F F -I- 9F F F F F^ p p2 p -p2 

J- li-L jj^ kk T iji- jk^ ki -L jj-L -L jji i fcfci jj 



where we have dropped the argument p^ of the Tijijp') appearing on right-hand side for ease 
of notation. 

The complex pole of each propagator is determined as the solution of 

-M- -"^•+Sf(K) = 0. (75) 

Writing the complex pole as 

M'^ = M'^ - iMT, (76) 

where M is the mass of the particle and F its width, and expanding up to first order in F 
around yields the following equation for Mf, 

, , , ImSf(Mf) (lmEf)'(M2) 

- + Re (Mf ) + , = 0- (77) 

l+(ReSfj {Mi) 

As before, in Eq. (77) the short-hand notation f'{p^) = d fip^)/ [dp^) has been used, and 
Mi denotes the loop-corrected mass, while mj is the lowest-order mass {i = h, H, A). 

While the Higgs-boson masses can in principle directly be determined from Eq. (77) 
by means of an iterative procedure (since appears as argument of the self-energies 
in Eq. (77)), it is often more convenient to determine the mass eigenvalues from a diag- 
onalization of the mass matrix in Eq. (71). In our numerical analysis (and in the code 
FeynHiggs 2 . 5) we perform a numerical diagonalization of Eq. (71) using an iterative Jacobi- 
type algorithm [80] . The mass eigenvalues Mj are then determined as the zeros of the function 
l^fip^) —p"^, where /i^(p^) is the ith eigenvalue of the mass matrix in Eq. (71) evaluated at p"^. 
Insertion of the resulting eigenvalues Mj into Eq. (77) verifies (to 0(r)) that each eigenvalue 
indeed corresponds to the appropriate (complex pole) solution of the propagator. We define 
the loop-corrected mass eigenvalues according to 

Mh, < Mh, < Mh,. (78) 

In our determination of the Higgs-boson masses we take into account all imaginary parts 
of the Higgs-boson self-energies (besides the term with imaginary parts appearing explicitly 
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in Eq. (77), there are also products of imaginary parts in ReE^f (M^^)). The effects of the 
imaginary parts of the Higgs-boson self-energies on Higgs phenomenology can be especially 
relevant if the masses are close to each other. This has been analyzed in Ref. [58] taking into 
account the mixing between the two heavy neutral Higgs bosons, where the complex mass 
matrix has been diagonalized with a complex mixing angle, resulting in a non-unitary mixing 
matrix. The effects of imaginary parts of the Higgs-boson self-energies on physical processes 
with s-channel resonating Higgs bosons are discussed in Refs. [58-60]. In Ref. [58] only the 
one-loop corrections from the t/t sector have been taken into account for the H-A mixing, 
analyzing the effects on resonant Higgs production at a photon collider. In Ref. [59] the full 
one-loop imaginary parts of the self-energies have been evaluated for the mixing of the three 
neutral MSSM Higgs bosons. The effects have been analyzed for resonant Higgs production 
at the LHC, the ILC and a photon collider (however, the corresponding effects on the Higgs- 
boson masses have been neglected). In Ref. [60] the t/h one-loop contributions (neglecting 
the t/h corrections) on the H-A mixing for resonant Higgs production at a muon collider 
have been discussed. Our calculation incorporates for the first time the complete effects 
arising from the imaginary parts of the one-loop self-energies in the neutral Higgs-boson 
propagator matrix, including their effects on the Higgs masses and the Higgs couplings in a 
consistent way. 

As described above, the solution for the Higgs-boson masses in the general case where 
the full momentum dependence and all imaginary parts of the Higgs-boson self-energies are 
taken into account is numerically quite involved. It is therefore of interest to consider also 
approximate methods for determining the Higgs-boson masses (often used in the literature) 
and to investigate in how far the results obtained in this way deviate from the full result. 
Instead of keeping the full momentum dependence in Eq. (71), the "p^ on-shell" approxi- 
mation consists of setting the arguments of the self-energies appearing in Eq. (71) to the 
tree-level masses according to (i, j = h, H, A) 



%{P') - E.,((m? + m|)/2). 



on-shell approximation: ' - 2 2\ /n\ C''^) 



In this way the Higgs-boson masses can simply be obtained as the eigenvalues of the 
(momentum-independent) matrix of Eq. (71). The "p^ on-shell" approximation has the 
benefit that it removes all residual dependencies on the field renormalization constants that 
cannot be avoided in an iterative procedure for determining the mass eigenvalues, see the 
discussion in Sect. 2.7. 

Instead of setting the momentum argument of the renormalized self-energies to the tree- 
level masses, in the "p^ = 0" approximation the momentum dependence of the self-energies 
is neglected completely = h,H,A), 

p2 = approximation: ^''^J^l] ^ ^''^^\ (80) 

In the "p^ = 0" approximation the masses are identified with the eigenvalues of Mn(0) 
(see Eq. (71)) instead of the true pole masses. This approximation is mainly useful for 
comparisons with effective-potential calculations and the determination of effective couplings 
(see below). The matrix Mn(0) is hermitian (and real and symmetric) by construction. 
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In order to study the impact of the imaginary parts of the Higgs-boson self-energies, it 
is useful to compare the full result with the "Im S = 0" approximation, which is defined by 
performing the replacement 

ImS = approximation: S(p^) ReS(p^) (81) 

for all Higgs-boson self-energies. Also this approximation results in an hermitian mass 
matrix. The comparison of our full result with the "p^ on-shell", the "p^ = 0" and the 
"ImE = 0" approximations will be discussed in Sect. 4. 



3.5 Amplitudes with external Higgs bosons 



In evaluating processes with external (on-shell) Higgs bosons beyond lowest order one has to 
account for the mixing between the Higgs bosons in order to ensure that the outgoing particle 
has the correct on-shell properties such that the S matrix is properly normalized. This gives 
rise to finite wave-function normalization factors.^ For the case of 2 x 2 mixing appearing in 
the rMSSM for the mixing between the two neutral CP-even Higgs bosons h and H, which 
is analogous to the mixing of the photon and Z boson in the Standard Model, the relevant 
wave function normalization factors are well-known, see e.g. Refs. [21,81]. An amplitude 
with an external Higgs boson, i, receives the corrections = h,H, no summation over 
hj) 



r. + 



7 r 



^ j) 



where the Tij denote the one-particle irreducible Higgs vertices, and 

l + ReS^(p=^) - Re 



— rrij + E 



'3J 



^p2 



J 



--Mf 



(82) 



(83) 



4) 



Mf - m) + S,,(M2) 

As before denotes the tree-level mass, while Mj is the loop-corrected mass. 

In the case of the cMSSM, the formulas above need to be extended to the case of 3 x 3 
mixing. This can easily be achieved using the results of Sect. 3.4. A vertex with an external 
Higgs boson, z, has the form (with /c all different, A; = h,H,A, and no summation 
over indices) 



Ti + 



ZijVj 



+ ZikVk + 



^5) 



where the ellipsis represents contributions from the mixing with the Goldstone boson and 
the Z boson, as discussed above. The finite Z factors are given by 

^ (86) 



1 + 



ReSlf)'(Mf) 



^The introduction of these factors can in principle be avoided by using a renormalization scheme where 
all involved particles obey on-shell conditions from the start, but it is often more convenient to work in a 
different scheme like the DR scheme for the field renormalizations described in Sect. 2.6. 
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M? 



ml 



Mf - ml 



where the propagators Ajj(p^), Ajj(p^) have been given in Eqs. (72) and (74), respectively. 
Using Eq. (85) with the Z factors specified in Eqs. (86), (87) and adding to this expression 
the mixing contributions of the Higgs bosons with the Goldstone bosons and the gauge 
bosons (see the discussion above) yields the correct normalization of the outgoing Higgs 
bosons in the S matrix. 

For later convenience we define a matrix Zn based on the wave function normalization 
factors. Its elements are given by (with Za = 1, i,j = h, H, A, and no summation over i) 



(Z„)i 



Zi Zij 



Some care is necessary in order to correctly identify the elements (Zn)jj (given in terms of 
the h,H,A states) with the corresponding mass eigenstates hi,h2,hs. To find the correct 
assignment, besides using Eq. (77) as described above, for mass-degenerate cases we also 
compute the matrix Zn for all possible permutations of the Higgs bosons involved in the 



Cij\ 



Here Cij is the 



m 



mixing and choose the permutation which minimizes 
general non-unitary) mixing matrix resulting from diagonalizing the full mass matrix.^ This 
procedure results in the matrix Zn that is obtained from the matrix Zn by a re-ordering of 
its rows. A vertex with an external Higgs boson hi is then given by 



(Zn)jir/i + (Zn)i2^H + (Zn)i3rA + 



^9) 



where the ellipsis again represents contributions from the mixing with the Goldstone boson 
and the Z boson. 



3.6 Effective couplings 

In a general amplitude with internal Higgs bosons, the structure describing the Higgs part 
is given by 

E (90) 

where the r^j- are as above the one-particle irreducible Higgs vertices, and the propagators 
Aij are given in Eqs. (72) and (74). For phenomenological analyses it is often convenient 
to use approximations of improved-Born type with effective couplings incorporating leading 
higher-order effects. There is no unique prescription how to define such effective coupling 

^The matrix C depends of course on the external momentum where it is evaluated. Since the depen- 
dence on is not very pronounced and we need C only to distinguish mass-degenerate cases, we choose 
the C evaluated at — A/| since the mass ordering ensures that the second-lightest Higgs boson is always 
involved in the degeneracy. 
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terms. One possibility would be to consider the matrix Zn, defined through Eqs. (88)-(89), 
as mixing matrix. The elements of the matrix Zn, however, are in general complex, so that 
Zn is a non-unitary matrix. Therefore it cannot be interpreted as a rotation matrix. If one 
wants to introduce effective couplings by means of a (unitary) rotation matrix, it is necessary 
to make further approximations. 

A possible choice leading to a unitary rotation matrix is the "p^ = 0" approximation, 
which is used in the effective potential approach. As before, we first consider the case of 
2x2 mixing relevant for the rMSSM. In the "p^ = 0" approximation defined in Eq. (80) the 
momentum dependence in the renormalized self-energies is neg lected, S(p2) ^ S(0), so that 
the derivative in Eq. (83) acts only on the term in the propagator factor. In this limit Zi 
simplifies to [50,82] 

= approximation, 2x2 mixing: Zi = — . (91) 

1 + ^ 

For the mixing between the neutral CP-even Higgs bosons h, H this yields Zh = Zh = 
cos^ Aa. This corresponds to an effective coupling approximation where the tree-level mixing 
angle a appearing in the couplings of the neutral CP-even Higgs bosons is replaced by 
«efr = a + Aa [50,82]. 

It is easy to verify that for the 3x3 mixing case Eq. (86) in the "p^ = 0" approximation 
simplifies to 

= approximation, 3x3 mixing: Zi = — , (92) 

1 + + Z^f, 

as a direct generalization of Eq. (91). 

The matrix Zn defined through Eqs. (88)-(89) goes over into a unitary rotation matrix 
Rn in this approximation, 

I Ru Rl2 Rl3\ 

(93) 

The matrix Rn diagonalizes the matrix Mn(0) arising from Eq. (71) in the "p^ = 0" ap- 
proximation. Rn can therefore be used to connect the mass eigenstates /ii,/i2,/i3 with the 
original states h, H, A, 

RnMn(0)Rt= Ml^^,^^ . (94) 
V Ml^^J 

We will discuss in this paper also the possibility of defining the effective couplings in the 
"p^ on-shell" approximation. The unitary matrix Un is then defined such that it diagonalizes 
the matrix Re (Mn(p^ on-shell)) arising from Eq. (71) in the "p^ on-shell" approximation and 
restricting to the real part of the matrix. This yields 

on-shell approx., 3x3 mixing: | ^2 | = Un 

on— shell 







Ru 


R12 


R 


= approximation, 3x3 mixing: 


Zn - 


Rn, Rn = R2I 


R22 


R. 






\R31 


R32 


R: 

















p2=0 












f/12 


f/l3 




|,Un^| 


f/21 


f/22 


f/23 






\U31 


f/32 


f/33 
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fMl 2 

/ n-i,P on— shell 

Un Re (M,(p2 on-shell)) Uj; = M^^^, | . (95) 

V M2^_^2 „,_3hell. 

The elements of Un, which can be chosen to be real, can be used to quantify the extent 
of CP- violation. (The same applies to Rn, which is real by construction.) For example, 
can be understood as the CP-odd part in hi, while the combination Ufi + f/fg corresponds 
to the CP-even part. The unitarity of Un ensures that both parts add up to 1. 

The elements of Un (or Rn) can be interpreted as effective couplings of Higgs bosons, 
which take into account leading higher-order contributions. As an example, we discuss here 
the effective couplings of the neutral MSSM Higgs bosons to SM gauge bosons and fermions. 

Beyond the lowest order in the cMSSM all three neutral Higgs bosons have a CP-even 
component, so that all three Higgs bosons have non- vanishing couplings to two gauge bosons, 
VV = ZZ, W~^W~ . The couplings normalized to the SM values are given by 

gh.vv = Uii sin(/5 - a) + Ui2 cos(/3 - a). (96) 

The coupling of two Higgs bosons to a Z boson, normalized to the SM value, is given by 

9h,h,z = Uis (Uji cos{P -a)- Uj2 sin(/? - a)) 

- Ujs {Uii cos(/5 -a)- Ui2 sin(/? - a)) . (97) 

The Bose symmetry forbidding any anti-symmetric derivative coupling of a vector particle 
to two identical real scalar fields is respected, QhihiV = 0. 

Concerning the decay into light SM fermions, we will compare in Sect. 4 below the full 
result based on the wave function normalization factors with the effective coupling approx- 
imation. In the latter approximation, the decay width of hi can be obtained from the SM 
decay width of the Higgs boson by multiplying it with 



where 



atuu = (Uii cos a + Ui2 sin a)/sp, = Ui3 c^/ (99) 

atdd = {-Uii sin a + Ui2Cosa)/cp, gl^^ = Ui3 sp/cp (100) 

for up- and down-type quarks, respectively. 

The results obtained by using effective couplings for simplified calculations of cross sec- 
tions or decay widths at fixed-order perturbation theory are inherently less precise than those 
from a full diagrammatic calculation at the same order. If effective couplings are employed, 
their limitations should be kept in mind. It will be shown below that for not too large val- 
ues of Mh± effective couplings evaluated with Un give results closer to the full calculation 
of Eq. (89) for the propagator corrections on external lines than those evaluated with Rn. 
On the other hand, it can be shown analytically that the effective couplings of the lightest 
Higgs boson evaluated with Un do not decouple to the SM limit for Mh± oo. Decoupling 
can only be achieved employing either the full calculation of Eq. (89) or effective couplings 
evaluated with Rn. 
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4 Numerical analysis 



Our results obtained in this paper extend the known results in the literature in various 
ways. The results for the Higgs-boson masses and couplings in the cMSSM available so 
far have been restricted to evaluations in the EP approach [53,54] (at one-loop, neglecting 
the momentum dependent effects) and to the RG improved one- loop EP method [55,56]. 
In Refs. [53,55,56] only corrections from the (s)fermion sector and the gaugino sector have 
been taken into account, and various non-logarithmic terms and momentum-dependent cor- 
rections have been neglected. A calculation taking into account also contributions from the 
gauge-boson and Higgs sector has been performed in Ref. [54], however (besides neglecting 
momentum dependent effects) using the parameter (see Eq. (16)) as input. Within the 
FD approach so far only the leading one-loop corrections had been evaluated, using the 
on-shell renormalization scheme [57]. Effects of imaginary parts of the one-loop contribu- 
tions to Higgs masses and couplings have mostly been neglected in the above results. Some 
effects induced by products of imaginary parts have been considered in Refs. [58-60], see the 
discussion in Sect. 3.4. 

Our results are based on the complete one-loop results in the cMSSM, taking into account 
the full dependence on the complex phases, the other MSSM parameters, and the external 
momentum. They involve a consistent treatment of all imaginary parts appearing in one-loop 
Higgs-boson self-energies that contribute to the Higgs-boson masses and the wave-function 
normalization factors of external Higgs bosons. Our one-loop results are supplemented by the 
dominant two-loop corrections in the FD approach, as described in Sect. 3.3. The higher- 
order corrected Higgs-boson sector has been evaluated with the help of the Fortran code 
FeynHiggs2.5 [23,44-46], see Sect. 5 below. 

4.1 Parameters 

In the context of a detailed phenomenological analysis of the cMSSM parameter space the 
existing constraints on CP- violating parameters from experimental bounds [83,84] are of in- 
terest. The complex phases appearing in the cMSSM are experimentally constrained by their 
contribution to electric dipole moments of heavy quarks [85] , of the electron and the neutron 
(see Refs. [86,87] and references therein), and of deuterium [88]. While SM contributions 
enter only at the three-loop level, due to its complex phases the cMSSM can contribute al- 
ready at one-loop order. Large phases in the first two generations of (s)fermions can only be 
accommodated if these generations are assumed to be very heavy [89] or large cancellations 
occur [90], see however the discussion in Ref. [91]. In the chargino and neutralino sector 
the three parameters Mi, M2 and yU. can be complex. However, there are only two physical 
complex phases since one of the two phases of Mi and M2 can be rotated away. One finds 
that in particular the phase yj^ is tightly constrained (in the convention where ipM2 = 0). 
The bounds on the phases of the third generation trilinear couplings, on the other hand, 
are much weaker. In order to illustrate the possible effects of complex phases we will show 
below results for ipM^ as well as ipM2 varied over the full parameter range. We will discuss 
the impact of the experimental constraints where appropriate. We treat the gluino mass 
parameter, which enters the observables discussed below only from two-loop order on, as 
real, M3 = m^. 
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Our numerical analysis has been performed for the following set of parameters (if not 
indicated differently): 

MsusY = 500 GeV, \At\ = \Ak\ = \Ar\ = 1000 GeV, 

= 1000 GeV, IM2I = 500 GeV, |Mi| = 250 GeV, m-g = 500 GeV, 
Mh± = 150 GeV, tan/3 = 5, 15, /i^^- = = 171.4 GeV [92]. (101) 

We do not consider higher values of tan (3, which in general enhance the SUSY contributions 
to the electric dipole moments. 

In order to evaluate the possible size of CP-violating effects in the Higgs sector in a 
conservative way we have chosen a relatively low value of Mh±. Parts of the investigated 
parameter regions are challenged by the Higgs search performed at LEP [2,3], depending in 
particular on the parameters of the t sector. It should be noted, however, that within the 
cMSSM the limits from the Higgs search are in general weaker than in the rMSSM, giving 
rise even to situations where no experimental lower bound on M^^ can be established at 
all [3,93,94]. 

Our calculation at the one-loop level is completely general, containing all complex phases. 
Concerning the numerical analysis, as explained above, we restrict ourselves to low or mod- 
erate values of tan/?. Therefore the effects arising from the h/h sector stay relatively small. 
Consequently we do not study the effects of complex phases from this sector, but focus on 
the phases of Xt, At, and of the gaugino mass parameters Mi and M2. 

4.2 Predictions for the mass and couplings of the lightest Higgs 
boson 

We begin with the predictions for the mass and couplings of the lightest neutral Higgs boson 
of the cMSSM, which are of particular interest in view of the existing experimental bounds 
and of the prospective high-precision measurement of the mass of a light MSSM Higgs boson 
at the LHC and the ILC. We first compare our full result with the approximations discussed 
in Sect. 3.4. Furthermore we investigate the effects of the phases of M2 and Mi. We then 
compare the predictions for the partial decay widths of all three neutral Higgs bosons to 
r leptons based on the wave function normalization factors defined in Sect. 3.5 with the 
effective coupling approximation. 

4.2.1 Comparison of the full result with approximations 

In Fig. 4 the cMSSM prediction for the mass of the lightest neutral Higgs boson, M/i^, is 
shown in the upper two plots, while the lower plot shows the coupling of hi to gauge bosons. 
The results are displayed as a function of the complex phase ipxt for l^tl = 700 GeV. 
The other parameters are chosen as specified in Eq. (101). Varying ^pxt leaves the t masses 
unchanged, so that the impact of the phase dependence is not masked by the purely kinematic 
effect of a change in the t masses. Our full result is compared with various approximations. 
In the upper left plot the full result for all sectors of the cMSSM is compared with the 
results taking into account only the effects of the /// sector (dot-dashed) and from the 
t/t + h/h sector (dashed). It should be noted that the asymmetry between the results for 
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at (fxt = and Lpxt = which amounts to about 8 GeV in this example, arises 
both from X^-dependent one-loop corrections (whereas the leading one-loop corrections 
in the limit Ma,Mh± ^ Mz depend only on the absolute value of X^, see e.g. Ref. [24]) 
and from two-loop contributions. For the parameters chosen in Fig. 4 there is a partial 
compensation between the phase variation at the one-loop and the two-loop contributions. 
The corrections beyond the /// loops, arising from the chargino/neutralino sector, the gauge- 
boson sector and the Higgs sector, can amount up to about 3 GeV. The /// contributions 
are clearly dominated by the contributions of the third generation quarks and squarks, 
with a maximum deviation of about 1 GeV for ipxt ~ ±7r. Effects at the sub-GeV level 
may be probed at the LHC and the ILC, where the anticipated precision for measuring 
the mass of a light Higgs boson is about 0.2 GeV (LHC) [8] and 0.05 GeV (ILC) [13-15]. 
For a discussion of theoretical uncertainties from unknown higher-order corrections and the 
parametric uncertainties induced by the experimental errors of the input parameters, see 
e.g. Refs. [41-43]. 

The upper right plot of Fig. 4 shows the difference between the full result and the "p^ on- 
shell", 'y = 0" and "ImS = 0" approximations defined in Eqs. (79)-(81). The "p^ = 
0" approximation yields results that differ from the full result by up to 1.5 GeV in this 
example, while the "p^ on-shell" approximation agrees with the full result to better than 
about 0.5 GeV. As explained above, the imaginary parts in the one-loop Higgs-boson self- 
energies arise only from kinematical thresholds, while the complex parameters enter only in 
combinations that are real. As a consequence, for the chosen set of SUSY parameters the 
self-energies entering the prediction for the lightest cMSSM Higgs mass develop imaginary 
parts only from loops involving SM fermions (except the top quark). The effects of neglecting 
the imaginary parts are therefore very small in this example, and the result in the "Im S = 0" 
approximation is indistinguishable in the plot from the full result. 

The coupling of the lightest cMSSM Higgs boson to gauge bosons normalized to the 
SM Higgs boson coupling, |(7hiyyp (obtained using the "p^ on-shell" approximation, see 
Eq. (96)), is shown in the lower plot of Fig. 4 for the same set of parameters. This cou- 
pling governs the Higgs production cross section in the Higgs-strahlung channel at LEP, 
the Tevatron and the ILC as well as the weak-boson fusion cross section at the LHC. The 
full result incorporating the contributions from all sectors of the MSSM (full line) differs 
from the result based on the /// sector only (dot-dashed) by up to 5 (10)% in the case 
of tan/? = 5 (15). The fact that the contribution from the t/t + b/b sector yields a better 
approximation of the full result for Ifi'hiyyP than the contribution from the whole /// sector 
is due to an accidental cancellation of contributions from different MSSM sectors. 

4.2.2 Dependence on the gaugino phases 

We now analyze the dependence on the gaugino phases ipMi and fM2- In Fig- 5 the depen- 
dence of the lightest cMSSM Higgs-boson mass on (pM2 is shown. The difference AM/^^ : = 
M/j^(all sectors) — Mh^{f / f sector), which is dominated by the chargino/neutralino contri- 
butions, is evaluated for three different values of jMs], jMal = 200,1000,2000 GeV (sohd, 
dashed, dot-dashed line). The other parameters are chosen as in Eq. (101), and all other 
complex phases are set to zero. The result including the full momentum dependence is given 
by the upper set of curves, while the "p^ = 0" approximation is given by the lower set. In 
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Figure 4: M/j^ and Isfhiyy P are shown as a function of ^pxt for = 700 GeV, tan /5 = 5 
and the other parameters as given in Eq. (101). 
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Figure 5: AM^-^ := M/t^(all sectors) — Mh^{f / f sector) is shown as a function of ipM-z for 
the full result and the "p^ = 0" approximation. The left plot shows the result for tan /9 = 5, 
while in the right plot tan/5 = 15. jMal is chosen as 200, 1000, 2000 GeV. 
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Figure 6: AMh^ := M/i^(all sectors) — Mh^{f / f sector) is shown as a function of v?Mi for 
the full result and the "p^ = 0" approximation. The left plot shows the result for tan /3 = 5, 
while in the right plot tan/3 = 15. |Mi| is chosen as 200, 1000, 2000 GeV. 
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the left plot we have chosen tan /? = 5, in the right one tan [3 = 15. For the lower tan (3 value 
the effects from the non-sfermionic sector are about 2-3 GeV if the full dependence on the 
external momentum is taken into account, and about 4 GeV in the "p^ = 0" approximation 
(which is used in the effective potential approach). The effect arising from varying the gaug- 
ino phase (pM2 itself is of (9(1 GeV). Both the overall effect from the non-sfermionic sector 
and the effect from varying (pM-z become smaller for larger tan/3 values (right plot). The 
effects are largest for IM2I =0{1 TeV), i.e. for IM2I being of the same order as In this 
case the gaugino-higgsino mixing in the chargino and neutralino sector, and correspondingly 
the couplings of the charginos and neutralinos to the Higgs sector, is maximized. The effects 
shown in Fig. 5 arising from varying ipM2 should be interpreted as an upper bound on the 
possible impact of the phase dependence. The possible effects from the gaugino phases will 
be reduced if the existing experimental constraints on these phases are taken into account, 
see the discussion above. 

We now turn to the effects from varying (pMi as shown in Fig. 6. The parameters are as 
in Fig. 5, but with M2 = 500 GeV and \Mi\ = 200, 1000,2000 GeV (solid, dashed, dotted 
line). The size of the effects from the non-sfermion sector is the same as in Fig. 5. However, 
the dependence on ipMi is much smaller, being of 0(100 MeV). 

4.2.3 Decay widths of the neutral Higgs bosons 

In this section we compare the predictions for the partial decay widths of all three neutral 
Higgs bosons to r leptons based on the wave function normalization factors as given in 
Eqs. (88), (89) with the effective coupling approximation (using the "p^ on-shell" approx- 
imation) according to Eqs. (94), (95), and with the "p^ = 0" approximation as given in 
Eqs. (93), (94). In Fig. 7 we show 

r,,, := T{h, ^ r+r") and T{h, ^ r+r^)R, T{h, ^ r+r-)u (102) 

for i = 1,2,3 (upper, middle, lower row), where F refers to the full result based on the 
wave function normalization factors, and Fu, Fr correspond to the effective coupling ap- 
proximation evaluated with Un and Rn, respectively. Since we are only interested here in 
the comparison of the wave function normalization factors with the effective coupling ap- 
proximations, we omit the contributions arising from the mixing of the physical Higgs states 
with the Goldstone boson and the Z boson and we also do not take into account irreducible 
vertex corrections to the hiT^T~ vertices. The results are shown for Mu± = 150 GeV and 
\Xt\ = 700 GeV as a function of (fixt, where the other parameters are chosen according 
to Eq. (101) with tan/3 = 5 (left) and tan/3 = 15 (right). As a general feature it can be 
observed that Fu is closer to the full result F than Fr with only few exceptions (due to 
accidental numerical cancellations)^. This shows that the effective coupling defined through 
Un, Eq. (95), gives a somewhat better numerical description than the one defined through 
Rn, Eq. (94), as used in the effective potential approach. For tan/3 = 5 the deviations be- 
tween the "p^ = 0" approximation and the full result are mostly at or below the 5% level, 
where the largest effects in general occur in the decay width of the lightest Higgs boson. For 
tan /3 = 15 the absolute and relative deviations between the effective coupling approximation 

^For large values of M^± due to the non-decoupling effects in Fu, see the discussion in Sect. 3.6, Fr 
would give results closer to the full evaluation. 
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and the full result can be significantly larger in the case of the lightest Higgs boson. For 
the decay widths of hi the full result can differ from the "p^ = 0" approximation by more 
than 10%. In particular, in the region where T{hi t^t~) has a minimum {(fxt ~ i^) the 
relative deviation between the full result and reaches more than 25%. Also in this case 
the deviation between the full result and Fu, based on the "p^ on-shell" approximation, is 
much smaller. The deviations for h2 and are again at the level of 5%. For larger values of 
\Xt\ even larger differences between the effective coupling approximation and the full result 
can be found. 

4.3 Mass difference and mixing of the heavy neutral Higgs bosons 

We now turn to the predictions for the masses and the mixing of the heavy neutral Higgs 
bosons of the cMSSM. The discovery of heavy Higgs bosons (in addition to a light one) 
would clearly establish an enlarged Higgs sector as compared to the SM. In the cMSSM the 
two heavy neutral Higgs bosons /i2 and are in general relatively close in mass, so that the 
mixing induced by the CP-violating phases can give rise to resonance-type effects. We first 
analyze the mass difference of the two heavy neutral Higgs bosons, AM32 := M/13 — M/ij, 
in scenarios where the Higgs-boson self-energies can be enhanced by threshold effects. We 
then investigate the phase dependence of AM32 and discuss in how far this observable can 
be employed for distinguishing the cMSSM from the rMSSM. Finally we perform a detailed 
analysis of the mixing of h2 and that is induced by the presence of complex phases. 

4.3.1 Threshold effects for heavy Higgs bosons 

We first analyze the effects of thresholds appearing in the Higgs-boson self-energy diagrams 
(e.g. for rriA = m^^ +^{3' sixth diagram in Fig. 1) on the masses of the heavy neutral 

Higgs bosons. In the first two lines of Fig. 8 the mass difference AM32 := — M^^ is 
shown as a function of Mh± for tan/3 = 5 (left) and tan/3 = 15 (right) for two different 
values of the phase of At, fAt = (upper row) and ipAt = 7r/2 (middle row). The other 
parameters are chosen as in Eq. (101). We compare the full result (solid lines) with the 
"p^ on-shell" (dot-dashed), "p^ = 0" (dashed) and "ImS = 0" (dotted) approximations 
defined in Eqs. (79)-(81). It can be seen for the full result that the threshold effects may 
lead to a significant enhancement of the mass splitting between the states /12 and h^, so 
that mass differences in excess of AM32 = 10 GeV can occur even for Mh± values in the 
decoupling region where Mh± ^ Mz- This behaviour is not reproduced in the "p^ = 0" 
approximation (which is used in the effective potential approach). On the other hand, it turns 
out that the "p^ on-shell" approximation, see Eq. (79), gives a rather good approximation 
to the full result. The remaining deviations stay below the level of 1 GeV. It should be 
noted in this context that the sharp peaks displayed in Fig. 8 would get smoothened if the 
effects of finite widths of the internal particles in the Higgs-boson self-energies were taken 
into account. A precise prediction directly at threshold would require a dedicated analysis 
that is beyond the scope of the present paper. 

We now turn to the effects of the imaginary parts of the Higgs-boson self-energies, i.e. 
the comparison of the full result with the "ImS = 0" approximation as defined in Eq. (81). 
While for the example of the lightest cMSSM Higgs boson, shown in Fig. 4, the result for 
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Figure 7: The decay widths T{hi t^t~), T{hi T^r~)u and T{hi t^t^)r. (see text) 
are shown for z = 1, 2, 3 (upper, middle, lower row) as a function of (pxt with \Xt \ = 700 GeV. 
In the left column tan/3 = 5, while for the plots in the right column tan/3 = 15. The other 
parameters are chosen according to Eq. (101). 
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Figure 8: The mass difference AM32 := M/13 — is shown as a function of Mh± (upper 
and middle rows) and as a function of ip^t and tan (3 (lower row) for the parameters given in 
Eq. (101). The upper row shows the results for ip^^ = with tan/3 = 5 (left) and tan/5 = 15 
(right). The middle row shows AM32 for if At = 7r/2 with tan/3 = 5 (left) and tan/3 = 15 
(right). The lower row shows AM32 for tan/3 = 5, 15 and Mh± = 1000 GeV as a function of 
(fAt (left) and for ipAt = vr, Mh± = 700, 1000 GeV as a function of tan/3 (right). 
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neglected imaginary parts was indistinguishable from the full result, for the mass difference 
of the heavy Higgs bosons a difference becomes visible in the two upper lines of Fig. 8 around 
the thresholds. The effect is shown in more detail in the lower line of Fig. 8. In the lower left 
plot we show AM32 as a function of ipAt for t^'ii /3 = 5, 15 and Mh± = 1000 GeV. In the right 
plot we display AM32 as a function of tan/3 for Mh± = 700, 1000 GeV and ipAt = t^- The 
other parameters are again those of Eq. (101). As one can see from the plot, the difference 
between the full result and the approximation with neglected imaginary parts can be as large 
as about 5 GeV. 

4.3.2 Phase dependence of AM32 

We now analyze the dependence of the mass difference AM32 := M^.^ — M^^ on the complex 
phases in more detail, in particular in view of the question whether the detection of a certain 
mass difference between the two heavy neutral Higgs bosons could be a direct indication 
of non-zero complex phases in the MSSM. In Fig. 9 we show AM32 in the KeXf-lmXt 
plane for tan/3 = 5 (left) and tan/3 = 15 (right) for Mh± = 150 GeV (upper row) and 
Mh± = 500 GeV (lower row). The other parameters are given in Eq. (101). 

The results in Fig. 9 show that the smallest mass differences between M^g and occur if 
Xt = At — ^/ tan (3 is real or has only a relatively small imaginary part. For Mh± = 150 GeV 
and low tan (3 this happens only around Xt ~ 1000 GeV, while for higher tan /3 three minima 
are reached for Xt ~ —1200,-200,800 GeV. The largest mass differences are realized for 
relatively large \Xt\. While for tan /? = 5 all possible mass differences occurring in the Re Xt- 
ImXf plane are also realized on the real axis, for tan/3 = 15 the largest mass differences can 
only be found for a non-zero imaginary part of Xt- 

The qualitative behaviour changes somewhat for Mh± = 500 GeV. Again the smallest 
mass differences between M^^ and M^^ occur for small imaginary parts of Xt. Two minima 
are found for tan /? = 5. On the other hand, for tan /3 = 15 only the region of the Re Xj-Im Xt 
parameter space around Xt = Q results in a small value of AM32. The rather symmetric 
shape of the plot for tan/3 = 15 and Mh± = 500 GeV around Xt = shows that in this 
case the dominant contribution to AM32 depends only on the absolute value \Xt\ ~ \At\. 
For tan /3 = 5, on the other hand, the minimum values of AM32 are reached for both 
lie At ^ and KeXt 7^ 0. Similarly to the case of Mh± = 150 GeV and low tan/3 we find 
also for Mh± = 500 GeV (both for low and high tan /3) that a large mass difference AM32 
does not necessarily require a non-zero complex phase of Xt. Indeed, for large Mh± all mass 
differences realized for a parameter point in the complex Xt plane are also realized on its real 
axis. This means that the determination of the mass difference AM32 alone will in general 
not be sufficient to obtain direct information about the size of the complex phases. On the 
other hand, the interpretation of an observed mass difference in terms of the underlying 
SUSY parameters will be different in the rMSSM and the cMSSM. Valuable information 
for determining the parameters of the cMSSM including their complex phases can therefore 
be obtained by combining the mass difference AM32 with a suitable set of observables that 
exhibit a non-trivial dependence on the complex phases. 

We have investigated also the effects of the complex phases of M2 and Mi on the mass dif- 
ference AM32. The effects stay below the 1 GeV level for most parts of the parameter space. 
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Figure 9: AM32 := M^.^ — M^^ is shown in the ReXj-ImX^ plane for tan/3 = 5 (left) and 
tan/3 = 15 (right) for Mh± = 150 GeV (upper row) and Mh± = 500 GeV (lower row). The 
other parameters are as given in Eq. (101). 



4.3.3 Mixing of the heavy neutral Higgs bosons 

We finally analyze the mixing of the two heavy neutral Higgs bosons. Mixing effects, es- 
pecially between the second and the third Higgs boson, can potentially be sizable in the 
cMSSM. In Ref. [54] it was argued that the phases of the gaugino mass parameters play an 
important role in this context. 

As explained above, the elements of Un (or Rn) can be used to quantify the extent of 
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CP-violation. For example, the combination + can be understood as the CP-even 
part in h^, while the component U^.^ corresponds to the CP-odd part. The unitarity of Un 
ensures that both parts add up to 1. 

We first discuss the CP-conserving case, where U^^ is either 1 or 0, depending on the mass 
ordering of Mh and Ma (the higher-order corrected masses of the heavy neutral CP eigen- 
states). In Fig. 10 we show the mass difference AM32 = M^.^ — M^^ = \Mh — Ma\ together 
with the effective couplings t/fg (based on the "p^ on-shell" approximation, see Eq. (95)) 
and i?|3 (obtained in the "p^ = 0" approximation, see Eq. (94)). The three quantities are 
given as a function of Xf, which in the CP-conserving case is a real parameter. We have set 
Mh± = 500 GeV, tan/5 = 5 (tan/3 = 15) in the left (right) plot, and the other parameters 
are chosen according to Eq. (101). The change from 1 to in the (33) element of the ro- 
tation matrix should obviously occur at the same value of Xf where the mass hierarchy of 
the states H and A is inverted, i.e. where AM32 = 0. This correlation between the masses 
and the effective couplings is not automatic, however, since the masses have been calculated 
using the full higher-order corrections, while as discussed in Sect. 3.6 the effective couplings 
can only be obtained using certain approximations. Fig. 10 shows that the behaviour of 
C/I3 with Xt is well matched to the one of AM32, i.e. the step in U^^ occurs very close to 
the Xt value where AM32 = 0. For i?33, on the other hand, the behaviour of the effective 
coupling significantly differs from the one of the higher-order corrected masses. This effect 
is particularly pronounced for small tan/3 as can be seen in the left plot, where the value of 
Xt for which AM32 ~ is reached differs by more than 100 GeV from the corresponding Xt 
value for which R^^ changes from to 1. For tan/3 = 15 (right plot) the deviation is smaller 
but still significant. Fig. 10 clearly shows that those contributions which are omitted if ef- 
fective couplings are constructed using the "p^ = 0" approximation (as done in the effective 
potential approach) can be numerically sizable and important for a physically well-behaved 
result. We find also in this case (for not too large values of Mh±) that a better numerical 
description is obtained with effective couplings defined through Un, Eq. (95). At large Mh± 
values both matrices are insufficient for a precise description. 

Now we focus on the mixing of the heavy neutral Higgs bosons in the presence of complex 
parameters. As a measure of the mixing between these two states we show in Fig. 11 U^^ 
in the Re X^-Im Xt plane. The choice for Mh± and tan (3 is the same as in Fig. 9, and the 
other parameters are specified in Eq. (101). We have checked that Uf^ is very close to zero 
(i.e. the lightest Higgs boson is nearly a pure CP-even state) and U^^ ^ 1 — The mixing 
varies strongly with ipxt for both values of tan/3 and low and high Mh±- In particular, 
for relatively large values of the variation of the phase tpxt (with \Xt\ kept fixed) can 
cause t/|3 (and consequently also f/|3) to take on any value in the range < t/|3 < 1. It 
can furthermore be seen in all four panels of Fig. 11 that a large mixing between the two 
heavy neutral Higgs bosons, corresponding to the parameter regions where t/|3 ~ 0.5, is a 
feature that can happen quite easily in the MSSM with complex parameters. Studying the 
properties of the heavy Higgs bosons is therefore of particular interest, since they could, at 
least in principle, give access to large CP-violating effects. 

The connection between AM32 (the mass difference between the two heavy Higgs bosons) 
and [/I3 (the mixing of h2 and h^, i.e. the CP composition of the two heavy Higgs bosons) 
can be analyzed by comparing Fig. 11 with Fig. 9 (the choice for Mh± and tan (3 is the same 
in both figures). The regions of the nodal points in Fig. 11, i.e. the points in which a change 
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Figure 10: The mass difference AM32 = \Mh — Ma\ and the effective couphngs U^^ and 
i?33 (based on the "p^ on-shell" and "p^ = 0" approximations, respectively) are shown as 
a function of Xt (chosen to be real) for Mh± = 500 GeV, tan /5 = 5 (left) and tan /3 = 15 
(right). The other parameters are chosen according to Eq. (101). 

in Xf causes the largest variation in U^.^, coincide with the regions where the mass difference 
AM32 is close to zero. This behaviour, which occurs for all Mh± and tan (3 values, is clearly a 
resonance- type effect: in the parameter regions where the masses of the Higgs states become 
degenerate, the mixing effects between the states are maximal. 

Finally we analyze the effect of the gaugino phases on the mixing of the heavy neutral 
Higgs bosons. The dependence of U^^ on the gaugino phases is depicted in Fig. 12 for 
tan/3 = 5 and Mh± = 500 GeV. We have chosen KeXt as KeXt = 655 GeV such that 
ImXt = corresponds to a resonance region where AM32 is close to zero, namely the right 
nodal point in the lower left plot of Fig. 11. For displaying the effects of the gaugino phases 
ipM2 ^-iid ifMi we have chosen in Fig. 12 three different values of the imaginary part of X-t, 
ImXf = 0,20,200 GeV (solid, dashed, dot-dashed lines). For the case of the nodal point 
where ImXt = a. very strong variation of U^-^ with both (pM2 (left plot) and ipMi (right 
plot) is observed, covering the whole allowed range < U^^ < 1. It should be noted that one 
encounters in this case a strong variation of U^^ with the gaugino phases even if ipM2, fMi 
restricted to relatively small values. However, for ImXt = 20 GeV, i.e. only slightly away 
from the nodal point, the dependence of U^^ on the gaugino phases is already much smaller. 
For ImXt = 200 GeV the variation of U^^ with ipM2 s-nd ipMi is numerically insignificant. It 
becomes apparent that the gaugino phases can have a strong impact on the mixing of the 
heavy Higgs bosons, but only directly on resonance. Outside of the resonance regions the 
effects of the gaugino phases are small. This is in contrast to Ref. [54], where it was claimed 
that a strong dependence of the Higgs mixing on (pM2 and (pMi were a general feature of the 
cMSSM. 



39 



ReXf/GeV 



ReXf/GeV 



ImXf 
GeV 




-1000 1000 -1000 1000 

Mh± = 150 GeV 




Figure 11: U^^ is shown in the KeXt-lmXt plane for tan/? = 5 (left) and tan/5 = 15 (right) 
for Mh± = 150 GeV (upper row) and Mh± = 500 GeV (lower row). The other parameters 
are as given in Eq. (101). 
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5 The code FeynHiggs 2 . 5: program features 



FeynHiggs 2 . 5 [23, 44-46] is a Fortran code for the evaluation of the masses, decays and 
production processes of Higgs bosons in the MSSM with real or complex parameters. In 
this section we give a short overview about its features. More detailed information about 
installation and use can be found in the Appendix. 

The calculation of the higher-order corrections is based on the Feynman-diagrammatic 
(FD) approach as outlined in the previous sections. At the one-loop level, it consists of a 
complete evaluation, including the full momentum and phase dependence, and as a further 
option the full 6x6 non-minimal flavor violation (NMFV) contributions for scalar quarks [95, 
96]. At the two-loop level all available corrections from the real MSSM have been included 
(see Refs. [41, 97, 98] for reviews). They are supplemented by the resummation of the leading 
effects from the (scalar) b sector including the full complex phase dependence [99]. 

The loop-corrected pole masses are determined as the real parts of the complex poles 
as described in Sect. 3.4. The imaginary parts of the Higgs-boson self-energies are fully 
taken into account. The masses are evaluated with two independent numerical algorithms. 
Deviations between the two methods indicate potential problems of numerical stability. In 
addition to the Higgs-boson masses, the program also provides results for the effective Higgs- 
boson couplings and the wave function normalization factors for external Higgs bosons as 
described in Sects. 3.5,3.6. 

Besides the computation of the Higgs-boson masses, effective couplings and wave function 
normalization factors, the program also evaluates an estimate for the theory uncertainties of 
these quantities due to unknown higher-order corrections. The total uncertainty is the sum 
of deviations from the central value, AX = Yl'i=i ~ with X = {Mh^ h^^h^ H±, Uij, Zij}, 
where Uij is defined in Eq. (95) and Zij in Eqs. (88)-(89). Alternatively instead of Uij also 
Rij, defined in Eq. (94), can be evaluated. The Xi are given by 

• Xi is obtained by varying the renormalization scale (entering via the DR renormaliza- 
tion) within ^rrit ^ fi ^ '^mt, 

• X2 is obtained by using m^°^'^ instead of the running rrit in the two-loop corrections, 

• X3 is obtained by using an unresummed bottom Yukawa coupling, y^, i.e. a yb including 
the leading 0{asab) corrections, but not resummed to all orders. 

Furthermore FeynHiggs 2 . 5 contains the evaluation of all relevant Higgs-boson decay 
widths and effective couplings (the latter are given in the conventions used in the MSSM 
model file of the program FeynArts [70]). In particular, the following quantities are calcu- 
lated: 

• the total width for the neutral and charged Higgs bosons, 

• the branching ratios and effective couplings of the three neutral Higgs bosons to 

— SM fermions (see also Refs. [50,97]), hi //, 

— SM gauge bosons (possibly off-shell), hi 77, ZZ*, WW*,gg, 

— gauge and Higgs bosons, hi Zhj, hi hjhk, 
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— scalar fermions, hi p f\ 

— gauginos, hi XkXj, hi xl'Xm, 

• the branching ratios and effective couphngs of the charged Higgs boson to 

— SM fermions, H~ ff, 

— a gauge and Higgs boson, H~ hiW~ , 

— scalar fermions, H~ f'' f, 

— gauginos, H~ ^ xlxt 

• the production cross sections of the neutral Higgs bosons at the Tevatron and the LHC 
in the approximation where the corresponding SM cross section is rescaled by the ratios 
of the corresponding partial widths in the MSSM and the SM or by the wave function 
normalization factors for external Higgs bosons as defined through Eqs. (88)-(89), see 
Ref. [100] for further details. 

For comparisons with the SM, the following quantities are also evaluated for SM Higgs bosons 
with the same mass as the three neutral MSSM Higgs bosons: 

• the total decay width, 

• the couplings and branching ratios of a SM Higgs boson to SM fermions, 

• the couplings and branching ratios of a SM Higgs boson to SM gauge bosons (possibly 
off-shell). 

• the production cross sections at the Tevatron and the LHC [100]. 

FeynHiggs 2 . 5 furthermore provides results for electroweak precision observables that give 
rise to constraints on the SUSY parameter space (see Ref. [42] and references therein): 

• the quantity Ap up to the two-loop level [101] that can be used to indicate disfavored 
scalar top and bottom mass combinations, 

• an evaluation of My/ and sin^ ^efr, where the SUSY contributions are treated in the Ap 
approximation (see e.g. Ref. [42]), taking into account at the one-loop level the effects 
of complex phases in the scalar top/bottom sector [102] as well as NMFV effects [95], 

• the anomalous magnetic moment of the muon, including a full one-loop calculation [103, 
104] as well as leading and subleading two-loop corrections [105,106] (see also Ref. [107]), 

• the evaluation of BR(6 — ^ 57) including NMFV effects [96] . 
Finally, FeynHiggs 2 . 5 possesses some further features: 

• Transformation of the input parameters from the DR to the on-shell scheme (for the 
scalar top and bottom parameters), including the full 0{as) and 0{at^b) corrections. 
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• Processing of SUSY Les Houches Accord (SLHA 2) data [108-110]. FeynHiggs2.5 
reads the output of a spectrum generator file and evaluates the Higgs boson masses, 
branching ratios etc. The results are written in the SLHA format to a new output file. 

• Predefined input files for the SPS benchmark scenarios [111] and the Les Houches 
benchmarks for Higgs boson searches at hadron colliders [112] are included. 

• Detailed information about all the features of FeynHiggs 2 . 5 are provided in man 
pages. 

6 Conclusions 

We have presented new results for the complete one-loop contributions to the masses and 
mixing effects in the Higgs-boson sector of the MSSM with complex parameters. They 
have been obtained in the Feynman-diagrammatic approach using a hybrid renormalization 
scheme where the masses are renormalized on-shell, while the DR scheme is applied for tan f3 
and the field renormalizations. A detailed description has been given of the renormalization 
procedure and of the determination of the masses as the real parts of the complex poles 
of the higher-order corrected Higgs propagator matrix. Besides the Higgs-boson masses we 
have also derived the wave function normalization factors needed for processes with external 
Higgs bosons. We have discussed different ways for defining effective Higgs-boson couplings 
that incorporate leading higher-order effects. As a result, we propose effective couplings 
based on the "p^ on-shell" approximation, where the Higgs-boson self-energies are evaluated 
at the tree-level masses. 

In our calculation of the Higgs-boson masses, couplings and wave function normalizations 
the full dependence on all relevant complex phases is taken into account. We incorporate for 
the first time the complete effects arising from the imaginary parts of the one-loop Higgs- 
boson self-energies in a consistent way. Our result for the complete one-loop contributions 
in the cMSSM is supplemented by all available two-loop corrections in the rMSSM and a 
resummation of the leading effects from the sbottom sector for complex parameters. 

In our numerical discussion we have first analyzed the impact of our results on the physics 
of the light Higgs boson, which is of interest in view of the current exclusion bounds and 
possible high-precision measurements of the properties of a light Higgs boson at the next 
generation of colliders. We first investigated the impact of the different MSSM sectors and 
the possible effects of the corresponding complex phases on the mass of the lightest neutral 
Higgs boson of the cMSSM Higgs, M/^^, and the coupling of the lightest Higgs to gauge 
bosons. The well-known result from the rMSSM that the bulk of the corrections to Mh^ 
arises from the fermion/sfermion sector is of course also reflected in the dependence on the 
associated complex phases. We find that the effects associated with the variation of ipxt ^^re in 
general numerically very important, leading to shifts in M^^ of up to 8 GeV in the examples 
that we have studied. The corrections beyond the fermion/sfermion loops, arising from the 
chargino/neutralino sector, the gauge-boson sector and the Higgs sector, can amount up to 
about 3 GeV. The dependence of Mh-^ on the gaugino phases ipM^ and (pMi is in general 
rather small and will be difficult to resolve even in the high-precision environment of the 
ILC (in particular in view of the existing experimental constraints on the gaugino phases). 
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We have furthermore compared our full result with various aproximations. We find sizable 
deviations of up to 1.5 GeV in M^^ between the full result and the "p^ = 0" approxima- 
tion, which is often used in the literature. We find that the "p^ on-shell" approximation is 
significantly closer to the full result with maximum deviations below 0.5 GeV. 

For the example of the partial decay widths of the three neutral Higgs bosons into r 
leptons we have compared the results based on the wave function normalization factors of 
the external Higgs bosons with effective-coupling approximations. The dependence of the 
partial widths on the phase ^pxt is very pronounced. We find that the result based on an 
effective coupling in the "p^ = 0" approximation (corresponding to the effective potential 
approach) deviates from the full result by typically up to 5%, with maximum deviations of 
more than 10%. The effective couplings in the "p^ on-shell" approximation show a better 
agreement with the full result for not too large values of Mh± ■ 

While over a large part of the cMSSM parameter space the lightest neutral Higgs boson 
is almost a pure CP-even state, large CP-violating effects may influence the masses and 
mixings of the two heavier neutral Higgs bosons of the cMSSM. We have analyzed the impact 
of complex phases on the mass difference between the two heavy Higgs bosons, AM32, and 
on their mixing properties. Our full result has been compared with various approximations. 
We find that the mass difference AM32 can significantly be enhanced by threshold effects, 
so that mass differences of more than 10 GeV are possible even in the decoupling region 
where Mh± 3> Mz- Since the threshold effects go beyond the "p^ = 0" approximation, large 
deviations between the full result and this approximation may occur. The "p^ on-shell" 
approximation, on the other hand, is close to the full result even in the threshold region. We 
have furthermore shown that effects of the imaginary parts of the Higgs-boson self-energies 
can be sizable in the threshold region. They can be as large as about 5 GeV in this region. 

While the mass difference AM32 is sensitive to the effects of the complex phases in 
the sfermion sector (the effects of the gaugino phases, on the other hand, are small), a 
determination of AM32 will in general not be sufficient to establish the existence of non-zero 
complex phases. We have shown that most values of AM32 that can be obtained in the 
complex KeXt-lmXt plane can also be realized on the real axis, i.e. for ImXt = 0. In 
order to extract information on the complex phases, the mass difference AM32 will have to 
be combined with a suitable set of observables that exhibit a non-trivial dependence on the 
complex phases. 

We find that a large (CP-violating) mixing between the two heavy Higgs states is possible 
over a significant part of the cMSSM parameter space. In the parameter regions where the 
mass difference AM32 becomes very small a resonance-type behaviour is possible. It gives 
rise to large variations in the mixing between the two Higgs bosons, i.e. a small change in 
the phase (fxt can have a dramatic effect on the mixing properties. Directly on resonance 
even the gaugino phases ipM2 and ipMi have a large impact on the Higgs mixing. In contrast 
to Ref . [54] , where it was claimed that a strong dependence of the Higgs mixing on (p M2 a-^d 
^PMl were a general feature of the cMSSM, we find that outside of the resonance regions the 
effects of the gaugino phases are small. For a reliable description of the resonance region it 
is crucial to correctly take into account the imaginary parts of the Higgs-boson self-energies. 

Our analysis has shown that effective couplings in the "p^ = 0" approximation, as used in 
the effective potential approach, can be insufficient for correctly matching the Higgs-mixing 
properties to the higher-order corrected Higgs-boson masses. The effective couplings based 
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on the "p^ on-shell" approximation that we have studied in this paper, on the other hand, 
have turned out to be well-suited for a numerical description of the Higgs-boson mixing for 
not too large values of Mh± . However, for large Mfj± the correct decoupling properties of the 
effective couplings of the lightest Higgs boson are achieved in the "p^ = 0" approximation, 
but not in the "p^ on-shell" approximation. 

Our results for the Higgs-boson masses, couplings and wave function normalization factors 
together with an estimate of the remaining theoretical uncertainties from unknown higher- 
order corrections are implemented into the public Fortran code FeynHiggs 2 . 5. The code 
also contains the evaluation of the Higgs-boson decays and the main Higgs-boson production 
channels at the Tevatron and the LHC, calculated using the full wave function normalization 
factors. Further quantities that are useful for deriving constraints on the SUSY parameter 
space are also evaluated, such as electroweak precision observables, the anomalous magnetic 
moment of the muon and (in the case of complex parameters) electric dipole moments. The 
code can be obtained from www.feynhiggs.de . 
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A The code FeynHiggs 2 . 5: installation and use 
A.l Installation and Use 

The installation process is straightforward and should take no more than a few minutes: 

• Download the latest version from www.feynhiggs.de and unpack the tar archive. 

• The package is built with ./configure and make. This creates the library libFH.a 
and the command-line frontend FeynHiggs. 

• To build also the Mathematica frontend MFeynHiggs, invoke make all. 

• make install installs the files into a platform-dependent directory tree, for example 
i586-linux/{bin, lib, include}. 

• Finally, remove the intermediate files with make clean. 
FeynHiggs 2 . 5 has four modes of operation, 

• Library Mode: Invoke the FeynHiggs routines from a Fortran or C/C++ program 
linked against the libFH.a library. 

• Command-line Mode: Process parameter files in native FeynHiggs or SLHA format at 
the shell prompt or in scripts with the standalone executable FeynHiggs. 

• WWW Mode: Interactively choose the parameters at the FeynHiggs User Control 
Center (FHUCC) and obtain the results on-line. 

• Mathematica Mode: Access the FeynHiggs routines in Mathematica via MathLink 
with MFeynHiggs. 

A. 1.1 Library Mode 

The core functionality of FeynHiggs 2.5 is implemented in a static Fortran 77 library 
libFH.a. All other interfaces are 'just' frontends to this library. 

In view of Fortran's lack of symbol scoping, all internal symbols have been prefixed 
to make symbol collisions very unlikely. Also, the library contains only subroutines, no 
functions, which simplifies the invocation. In Fortran, no include files are needed except for 
access to the coupling structure. In C/C++, a single include file CFeynHiggs.h must be 
included once for the prototypes. Detailed debugging output can be turned on at run time. 

The library provides the following functions: 

• FHSetFlags sets the flags for the calculation. 

• FHSetPara sets the input parameters directly, or 
FHSetSLHA sets the input parameters from SLHA data. 

• FHSetCKM sets the elements of the CKM matrix. 
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• FHSetNMFV sets the off-diagonal soft SUSY-breaking parameters in tlie scalar quark 
sector that induce NMFV effects. 

• FHSetDebug sets the debugging level. 

• FHGetPara retrieves (some of) the MSSM parameters calculated from the input pa- 
rameters, e.g. the sfermion masses. 

• FHHiggsCorr computes the corrected Higgs masses, effective couplings and wave func- 
tion normalization factors. 

• FHUncertainties estimates the uncertainties of the Higgs masses, effective couplings 
and wave function normalization factors. 

• FHCouplings computes the Higgs couplings and BRs. 

• FHHiggsProd calculates the Higgs-boson production cross-sections at the Tevatron and 
the LHC. 

• FHConstraints evaluates further electroweak precision observables. 

These functions are described in detail in their respective man pages in the FeynHiggs 
package. 

A. 1.2 Command-line Mode 

The FeynHiggs executable is a command-line frontend to the libFH.a library. It is invoked 
at the shell prompt as 

FeynHiggs inputfile [flags [scalef actor] ] 

where 

• inputfile is the name of a parameter file (see below). 

• flags is an (optional) string of integers giving the flag values, e.g. 40030211 (for 
details see the description of FHSetFlags in the man pages). If flags is not specified, 
40020211 is used. The fifth flag controls the evaluation of the effective couplings, 
where correspongs to the "p^ on-shell" and 4 to the "p^ = 0" approximation. 

• scalef actor is an optional factor multiplying the renormalization scale. It is used to 
determine the dependence on the renormalization scale, e.g. by varying scalef actor 
from . 5 to 2. 

FeynHiggs understands two kinds of parameter files: 

• Files in SUSY Les Houches Accord (SLHA) format. In this case FeynHiggs adds the 
Higgs masses, mixings and decay widths to the SLHA data structure and writes the 
latter to a file inputfile . f h. 

In fact, FeynHiggs tries to read each file in SLHA format first, and if that fails, falls 
back to its native format. 
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• Files in its native format, for example 



MT 


171.4 


MB 


4.7 


MW 


80.392 


MZ 


91.1875 


MSusy 


500 


MAO 


200 


Abs(M_2) 


200 


Abs (MUE) 


1000 


TB 


5 


Abs(Xt) 


1000 


Abs(M_3) 


800 



The syntax should be pretty self-explanatory. Complex quantities can be given either in 
terms of absolute value Abs(X) and phase Arg(X), or as real part Re(X) and imaginary 
part Im(X). Abbreviations, summarizing several parameters (such as MSusy) can be 
used, or detailed information about the various soft SUSY-breaking parameters can be 
given. 

Furthermore, it is possible to define loops over parameters, to scan parts of parameter 
space. For example, 

TB 5 25 5 

MAO 100 800 *2 

MSusy 500 1000 /3 

declares three loops: 

a) over tan (3 from 5 to 25 linearly in steps of 5 (i.e. 5, 10, 15, 20, 25), 

b) over Ma from 100 to 800 logarithmically in steps of 2 (i.e. 100, 200, 400, 800), 

c) over MsusY from 500 to 1000 linearly in 3 steps (i.e. 500, 750, 1000). 

The output is written in a human-readable form to the screen. Since this can be quite 
lengthy, a % is printed in front of all lines with 'non-essential' information, e.g. the 
details on couplings and decay widths. Thus to display only the 'essential' information, 
one just has to "grep" away the % lines, i.e. 

FeynHiggs inputfile flags I grep -v % 

The output can also be piped through the table filter to yield a machine-readable 
version appropriate for plotting etc. For example, 

FeynHiggs inputfile flags I table TB MhO > outputfile 

creates outputfile with two columns, tan/5 and M^. The syntax of the output file is 
given as screen output. 
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Debugging output is governed by the environment variable FHDEBUG which can be set to an 
integer from to 3 (for details see the description of FHSetDebug in the man pages). For 
example, 

setenv FHDEBUG 1 (in csh or tcsh) 

export FHDEBUG=1 (in sh or bash) 

sets debugging level 1. 
A.1.3 WWW Mode 

The FeynHiggs User Control Center (FHUCC) is a WWW interface to the command-line 
executable FeynHiggs. It provides a convenient way to play with parameters, but is of course 
not suited for large-scale parameter scans or extensive analyses. 
To use the FHUCC, point your favorite Web browser at 

http : / /www . f eynhiggs . de/f hucc 

adjust the parameters, and submit the form to see the results. At the end of the result page, 
the input file used for that FeynHiggs run is presented, too. 

A. 1.4 Mathematica Mode 

The MFeynHiggs executable provides access to the FeynHiggs functions from Mathematica 
via the MathLink protocol. This is particularly convenient both because FeynHiggs can be 
used interactively this way and because Mathematica's sophisticated numerical and graphical 
tools, e.g. FindMinimum, are available. 

After starting Mathematica, install the package with 

In[l]:= Install ["MFeynHiggs"] 

Out[l]= Link0bject[./i586-linux/bin/MFeynHiggs, 1, 1] 

which makes all FeynHiggs subroutines available as Mathematica functions. For details of 
their use, see the corresponding man pages. 
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